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ABSTRACT 


Full  wave  solutions  have  been  found  for  plane  wave  and  localized 
disturbances  in  the  ULF  (micropulsation)  and  ELF  frequency  ranges. 

Slow  and  fast  wave  disturbances  incident  at  5000  km.  altitude  at  the 
frequencies  100,  3.0,  1/3,  and  1/300  cps.  were  analyzed.  A  model 
ionosphere  developed  from  existing  rocket  and  satellite  data  for  conditions 
of  midday  and  midlatitude  at  sunspot  minimum  ..as  used  in  the  calculations. 
The  cold  plasma  dispersion  relation  was  written  conveniently  in  terms 
of  the  localized  disturbance  factor  and  the  ionospheric  conductivity 
tensor  elements.  The  localized  disturbances  were  found  to  produce  large 
vertical  electric  fields  in  the  neutral  atmosphere  for  the  incident  slow 
waves  in  the  micropulsation  range  and  the  incident  fast  wave  at  100  cps. 
Localization  causes  strong  mode  coupling  in  the  lower  D  region  and  the 
coupling  coefficients  are  larger  than  the  qurtic  roots  at  1/300  c/s. 

Charge  exchange  of  0+  and  O  results  in  an  enhanced  Pedersen  conductivity 
region  which  in  some  cases  acts  as  a  cavity  wail  for  wave  energy  storage. 
Charge  exchange  also  causes  some  absorption  at  Pul  frequencies  in  the 
range  150-250  km.  The  conductivity  of  the  Earth's  surface  is  included  and 
found  to  result  in  almost  total  reflection  at  ground  level  for  all  frequencies 
considered.  Incident  slow  mode  signals  at  3.0  cps.  undergo  polarization 
reversal  at  1000  km.  and  then  propagat  ..ostly  in  the  R  node  through 
the  ionosphere.  At  4000  km.  plasma  driif  velocities  ootween  0.1  and 
1 .0  km  /  sec.  are  sufficient  Hr  the  incident  sIoa  waves  to  excite  typical 
Pci  and  PcS  gn  and  observed  signals. 
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1.0  INTRODUCTION 


Variations  in  the  intensity  of  the  geomagnetic  field  as  observed 
at  the  Earth's  surface  are  well  known  to  exist  in  the  frequency  range 

-3 

1  cps.  down  to  around  10  cps.  These  oscillations  are  referred  to  as 
geomagnetic  micropulsations  and  have  been  given  a  morphological 
classification  by  Jacobs  et.  ai.  [1964].  It  is  now  generally  accepted 
that  micropulsations  as  observed  on  the  Earth's  surface  are  caused  by 
hydromagnetic  disturbances  generated  in  the  magnetosphere  and 
propagated  along  the  geomagnetic  field  lines  through  the  ionosphere 
down  to  ground  surface.  The  ultimate  aim  of  the  present  investigation 
is  to  determine  the  modification  of  these  waves  and  to  study  the  reflection 
and  magnetoionic  coupling  processes  as  the  disturbances  propagate 
down  to  ground  level.  A  full  wave  solution  is  obtained  utilizing  a 
numerical  technique  developed  by  Inoua  and  Horowitz  [1966a, b], 
Hydromagnetic  wave  propagation  has  been  studied  extensively 
ever  since  its  discovery  by  Alfven  [1942],  Theoretical  studies  led 
Dungey  [1954a,b]  to  suggest  that  some  micropulsations  are  of 
hydromagnetic  origin.  This  suggestion  has  since  led  to  a  number  of 
investigations  devoted  to  the  study  of  hydromagnetic  wave  propagation 
through  the  Earth's  ionosphere  [Akasofu,  1965;  Field,  1963;  Field  and 
Greifirgor,  1965,  I960;  Francis  and  Karplus,  1960;  Greifinger  and 
Grei  inger,  1965;  Jacobs  and  Wntanabe,  1962;  Karplus  et.  al.,  1962; 
Prlnco  and  9  s+ick,  1)64;  Prince  ut.  ai.,  1964], 


The  solutions  to  these  studies  have  been  effected  by  two  basic 
means,  i.e.,  analytical  and  numerical.  All  of  the  investigations 
mentioned  above  have  neglected  to  consider  magnetoionic  coupling, 
which  has  been  found  in  the  present  investigation  to  be  of  considerable 
importance  in  determining  the  final  wave  fields.  Analytical  solutions, 
such  as  those  obtained  by  Greifinger  and  Greifinger  [1965]  and  Field 
and  Greifinger  [1965,  1966],  have  been  effected  only  by  making  some 
drastic  simplifying  assumptions. 

For  high  and  medium  frequencies  the  wave  equation  is  often  solved 
using  the  WKB  method  and  the  approximations  of  geometric  optics .  For 
low  frequencies,  however,  the  concept  of  ray  optics  becomes  invalid 
and  partial  reflection  is  not  properly  treated.  A  more  accurate  solution 
to  the  wave  equation  must  then  be  sought  by  utilizing  known  functions 
or  numericai  methods;  this  is  referred  to  as  a  "full  wave  solution."  In 
general  full  wave  solutions  are  necessary  when  the  scale  size  of  the 
variations  in  the  ionospheric  parameters  becomes  less  than  the 
wavelength  of  the  electromagnetic  disturbance.  A  general  discus;  Ion 
of  the  full  wave  and  phase  integral  (i.e.,  WKB)  methods  and  thoi, 
advantages  and  disadvantages  may  be  found  in  an  excellent  revic.v 
article  by  Budden  [1969]. 

Various  full  wave  methods  have  teen  proposed  for  the  solution  of 
the  problem  of  electromagnetic  wave  propagation  in  a  horizontally 
stratified  magnetoionic  medium.  These  include  Budden  [1955],  Barron 
and  Budden  [1959],  Johler  and  Harper  [1962],  Pitteway  [1965],  anu 
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more  recently,  Walsh  [1967]  and  Altman  and  Cory  [I969a,b].  None 
of  these  methods,  with  the  exception  of  that  due  to  Altman  and  Cory, 
has  been  applied  to  the  micropulsation  frequency  range .  The  simple 
thin-film  optical  method  due  to  Altman  and  Cory  [1969a]  is  valid 
only  for  constant  wave  polarization,  when  magnetoionic  coupling  is 
zero.  In  general  this  condition  is  not  fulfilled  in  the  micropulsation 
frequency  range,  even  for  vertical  propagation.  The  generalized 
thin-film  optical  method  (ftltman  and  Cory,  [1969b])  allows  for  oblique 
incidence  and  intermode  coupling.  However,  the  wave  fields  do  not 
appear  to  have  been  calculated  with  this  method . 

In  the  present  investigation  three  different  frequencies  (3.0, 

0.333,  and  0.00333  cps.)  are  considered  in  the  micropulsation  frequency 
range.  In  addition  one  frequency  (100  cps.)  is  studied  in  the  ELF  range, 
which  is  of  interest  in  the  possible  effect  of  ion  cyclotron  resonance. 

For  each  frequency  a  localized  disturbance  (See  section  2.1)  is  considered. 
Also,  at  100  cps.  and  at  0.00333  cps.  plane  wave  cases  are  studied.  For 
each  case  considered  the  incident  slow  and  incident  fast  waves  are 
introduced  separately  at  the  top  of  the  ionosphere. 


2.0  THEORETICAL  CONSIDERATIONS 


2.1  Formulation  of  the  Wave  Equation 

The  assumption  is  made  that  an  electromagnetic  disturbance  is 
incident  at  5000  km.  altitude.  The  propagation  characteristics  of  this 
disturbance  ate  determined  bv  Max'.vell'c  equatir^s  and  the  constitutive 
relatione  for  the  medium.  The  dicturbat.ee  i.5  assumed  to  be  of 
sufficiently  small  amplitude  that  oouline-wiMes  of  tne  governing 
equations  may  be  neglected.  The  ionosphere  is  assumed  ‘o  be 
horizontally  stratified  in  the  presence  of  a  vertical  magnetic  field. 

The  disturbance  function  Fq  is  represented  by  a  double  Fourier 
series  over  the  spatial  coordinates  x  and  y. 

-L  v»\ 

where J£  and  m1  are  the  Fourier  numbers. 

The  Fourier  series  is  expanded  over  an  interval  equal  to  the  Earth's 
circumference  so  that  the  fundamental  wavelength  gives 

\  * 

where  =  Earth's  radius 

Maxwell's  equations  may  be  written  as 

vxe-  ik.# 

■7x#»-lk»R-E 
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\  a  free  space  wavelength 


tf  -  a.H 


where  ,  6^  ,  and  are  the  characteristic  impedance ,  dielectric 
constant,  and  magnetic  permeability  of  free  space  and  E  and  H  are 
the  electric  and  magnetic  wave  fields,  respectively.  The  quantity  W 
has  the  same  units  as  E  .  is  the  dielectric  tensor  of  the 
ionosphere  and  is  given  by 

-.1+1 
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where  X  is  the  unity  matrix  with  the  diagonal  elements  being  unity 
and  the  off-diagonal  elements  being  zero.  The  conductivity  tensor  (T 
is  determined  from  the  equations  of  motion  of  the  charged  particles 
(See  Part  II) . 

The  x  and  y  components  of  the  nabla  differential  operator  may  be 
written  as 

1 V- 


If  r  is  defined  as 
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then  Maxwell's  equations  in  component  form  may  be  written  as 


-TT  ^ 


=r^  =  <2/«' 


a 

-L-i-E-ie  ’r-W, -W' 
c  <  *  a  a 

i,  -m  ^ 

**>-»** 

$'  *  fV 

*'  -^ir 

These  equations  may  be  written  more  concisely  by  introducing  a 
itrix  notation  due  to  Clemmow  and  Heading  [1954]. 

e,-V 


rhere 


t  = 


*4 

•m; 

•5 


where  is  the  4  x  4  propagation  tensor  having  the  functional  form 

t»  r'tK'ijw) 

and  is  the  same  as  equation  00)  given  by  Inoue  and  Horowitz  [1966a] 
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with  K  being  replaced  by  . 

Although  the  propagation  tensor  T*  is  here  represented  as  a 
function  of  both  Hand  m,  a  coordinate  transformation  is  performed  so 
that  the  disturbance  is  propagating  only  in  the  x  and  z  directions . 

Then  the  Fourier  numbers  become 

Z  -  t*2  +  m'2)* 

m  =  o 

Also,  in  order  to  simplify  the  calculations  only  a  single  Fourier  component 
is  considered  so  that  the  disturbance  function  is  now  represented  by 
FD  ~  F  (>)  exp  ftltg*x)  e'1  Wt 

F  (s)  is  the  height  dependent  portion  of  the  disturbance  function  which 
is  determined  from  the  solution  to  the  wave  equation. 

The  above  differential  equation  represents  the  wave  equation 
reduced  to  a  set  of  four  linear,  homogeneous,  first  order  differential 
equations .  This  equation  was  Integrated  using  the  technique 
developed  by  Inoue  and  Horowits  [1966a  »b)  from  deep  in  the  Interior 
of  the  earth  out  to  SOOO  km .  altitude . 

In  order  to  carry  out  the  numerical  solutions  to  this  problem  a 
model  ionosphere  was  developed  Iran  existing  rocket  and  satellite 
data  for  conditions  of  mldiatitude  and  midday  at  sunspot  minimum. 

(See  Part  in  for  the  details  of  this  model.)  Hie  neutral  atmosphere 
fo  -  50  km.)  was  considered  to  be  free  space  with  zero  conductivity. 

The  surface  conductivity  of  the  earth  was  taken  from  data  published 
by  Green  and  Hoffman  [1967]  and  the  conductivity  of  the  Earth's 
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magnitude  of  cyclotron  frequency 
of  k  th  particle 


CO  =  angular  wave  frequency  of  the  disturbance 
The  quantities  n^,  m^,  z ^  represent  the  particle  density,  mass,  and 
charge  multiplicity  of  the  kth  particle  species.  In  this  work  z^  is 
assumed  to  be  unity,  e  is  the  magnitude  of  the  electrical  charge  of  an 
electron.  ^  is  +1  or  -1 ,  if  the  kth  particle  is  positively  or 
negatively  charged,  respectively.  BQ  is  the  Earth's  static  magnetic 
field . 

««X  -X 

The  curl  of  a  vector  A  is  57  x  A  and  may  be  written  as  an 
antisymmetric  tensor  times  the  vector  where  the  tensor  is  given  by 

°  ^ 

Ac  °  -$T 

Because  of  the  coordinate  transformation  which  mokes  Z  = 
and  m  =  o,  none  of  the  variables  have  any  y  dependence, 
derivatives  become 

h 

is  -*■ 


(£.2  +  m,2)4 
Tnen  the 
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where  q  is  the  Booker  quartic  root  (See  Budden,  1961,  p.  120). 


side  of  the  equation 
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=  0 


In  order  that  the  solution  to  this  system  of  equations  be 
nontrivial  the  determinant  of  the  3x3  matrix  must  be  zero.  The 
following  quartic  equation  representing  the  dispersion  relation  is 
then  obtained . 

Aq4  -  Bq2  +  C  =  0 

where 

Ar  5 

1  r  VlTfte  "]*“ 

rri  ~xr  1 


P,R,L  may  be  written  in  terms  of  C£  ,  0^  ,  and 
Hall,  and  Pedersen  conductivities,  respectively, 
section  4.1.) 


,  the  parallel , 
(See  Part  II, 
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Comparing  the  3  x  3  matrix  from  which  the  dispersion  relation 
is  derived  with  equation  (1-20)  given  by  Stix  [1962],  q  is  seen  to  be 
related  to  n,  the  total  index  of  refraction,  by 

M  /Cw,  9 
V 


where  Q  is  the  angle  between  the  wave  normal  and  the  magnetic 
field.  We  also  have  the  relation 

9  =  ^  ] 

Thus  we  see  that 

X  «  rw  Mw  s 


Now  0  can  be  eliminated  to  give  the  relationship  between  q,  n,  and  1. 

*  I  - 


<12  =  r2„2  - 
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Ifi=  o,  i.e.,  completely  vertical  propagation,  then 


Also,  in  this  case, 


t  =  TRL 

and  the  dispersion  relation  becomes 

o  9  2 

Substituting  q  =  r*  n  we  have 

Pn4  -  2PSn 2  +  PRL  =  O 

which  is  identical  to  the  dispersion  relation  given  by  Stix  for  0=0. 

2.3  Computer  Program  for  Full  Wave  Solution 

The  computer  program  for  the  full  wave  solution  consists  of  three 
separate  programs  (Part  I,  Part  II,  and  Part  III)  which  are  interfaced  by 
means  of  magnetic  tapes .  Quantities  calculated  in  Part  I  which  are 
necessary  for  the  computations  in  Part  II  are  stored  on  the  binary  magnetic 
tape  which  interfaces  Part  I  and  Part  II.  The  same  situation  is  true  for 
the  relationship  between  Part  II  and  Part  III. 

The  integration  procedure  used  in  this  program  is  based  on  the 
theoretical  work  of  Inoue  and  Horowitz  (1966a, b)  and  these  two  papers 
shall  hereinafter  be  referred  to  as  Papers  I  and  II,  respectively. 

2.31  Program  for  Part  I 

Part  I  of  the  full  wave  program  has  the  following  primary  functions: 
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TABLE  2.31 .  Some  important  variables  in  Part  I  and  their  meaning 


1 .  Performs  a  coordinate  transformation  so  that  the  propagation 
vector  lies  in  the  x  -  z  plane  at  an  angle  ©  to  the  z  axis. 

2.  Calculates  the  magneto- ionic  quanties  X,  Y,  Z. 

3 .  Uses  the  coupled  equations  of  motion  for  electrons  and 
ions  (See  Part  II,  chapter  1  of  this  report)  to  calculate  the  ionospheric 
conductivity  tensor. 

4 .  Calculates  the  dielectric  and  propagation  tensors . 

5.  The  altitude  (or  depth),  the  Fourier  number and  the 
dielectric  and  propagation  tensors  are  stored  on  magnetic  tape  for 
later  use  in  the  full  wave  program . 

A  few  of  the  important  variables  used  in  Part  I  are  given  in  Table 
2.31  along  with  their  interpretation.  The  wave  data,  station  data,  and 
medium  data  are  read  into  Part  I  on  cards . 

A  brief  description  of  the  various  subroutines  used  in  Part  I  is 
given  below . 

TABLE  2.32 


Subroutine 

DTRM 


DETERM 


MTMLTP 


MATINV 


Description 

Finds  the  determinant  of  a  3  x  3  matrix  having  complex 
elements . 

Computes  the  determinant  of  a  4  x  4  matrix  having 
complex  elements  using  the  method  of  pivot  condensation . 
Finds  the  product  of  two  3x3  matrices  each  having 
complex  elements . 

Determines  the  Inverse  of  a  3  x  3  matrix  having 
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isfjis 


complex  elements . 


DUMP  A  dump  is  forced  if  the  program  finds  itself  in  the 

wrong  location. 

FTRAP  A  University  of  Pittsburgh  subroutine  which  sets  to 

zero  any  number  which  becomes  too  small  for  the 
computer. 

Magnetic  tapes  used  in  the  program  are 

INP  BCD  scratch  tape  which  reads  in  data  cards 

IOUT  BCD  scratch  tape  which  prints  out  desired  information 

IW1  I  Binary  Tapes  used  for  manipulating  and  storing 

IW2  quantities  which  are  to  be  used  in  Part  II  and/or  Part  III . 

2.32  Program  for  Part  II 

Part  II  of  the  full  wave  program  has  the  following  primary  functions: 

1 .  A  binary  search  is  initiated  if  the  invariant  test  fails  (Paper  I, 
p.  433) .  This  test  will  fail  when  a  singularity  occurs  within  a  given  slab. 
The  purpose  of  the  binary  search  is  to  try  to  squeeze  the  singularity  into 
as  small  a  slab  as  possible. 

2 .  The  slab  width  is  adjusted  so  that  the  difference  in  the  amplitude 
changes  of  the  slow  and  fast  modes  within  a  slab  is  not  too  great. 

3 .  The  slab  is  determined  to  be  Isotropic  or  anisotropic  and  the 
Booker  quartic  roots  and  eigenvectors  are  calculated  accordingly. 

4.  The  coupling  coefficients  $ee  Paper  I,  p.  431)  are  calculated 
and  stored  on  tape  IW4,  which  is  used  as  a  permanent  storage  tape  for 
all  quantities  which  are  later  retrieved  and  plotted  automatically. 
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The  following  variables  are  read  into  Part  II  from  a  data  card. 

TABLE  2.33 

Variable  Description 

HMAX  Maximum  data  height  (km.).  For  purposes  of  the  work 

here,  HMAX=5000  km. 

HMIN  Minimum  data  height  (km.).  In  order  to  take  into 

account  the  conducting  properties  of  the  earth  this 
value  was  set  to  HMIN=-400  km. 

EFSINV  Numerical  value  assigned  to  the  invariant  test. 

WLNGTH  jqj  _  Re  -  Earth's  radius 

LIMBSR  The  limit  on  the  number  of  binary  searches  to  be 

preformed  for  a  given  slab. 

A  brief  description  of  the  various  subroutines  used  in  Part  II,  but 

not  in  Part  I,  is  given  below. 

TABLE  2.34 

Subroutine  Description 

COEFIC  Computes  the  coefficients  of  the  Booker  quartlc 

equation  from  the  propagation  tensor. 

POLAR  Changes  the  cartesian  coordinates  x  and  y  to 

polar  coordinates  r  and  $. 

LEIGEN  Calculates  the  normalised  matrix.  (See  Paper  I, 

p.  434). 

CMEXP  Given  the  real  and  imagiuary  parts  xr  and  rk, 

Xr+iXj 

respectively,  of  the  variable  x,  calculates  e. 


1-17 


BOOKER 


Calculates  the  roots  of  the  Booker  quartic  equation. 


Finding  a  root  requires  making  a  good  initial  guess 
for  the  root,  as  is  done  in  the  following  method.  The 
quartic  equation  is  differentiated  three  times  and  the 
zero  crossing  of  the  3rd 
derivative  is  used  to  find 
the  maximum  or  minimum 
of  the  2nd  derivative .  A 
Taylor  series  expansion  is 
then  made  about  the  maximum 


or  minimum  to  find  the  roots  of  the  2nd  derivative. 
Newton's  method  is  used  to  improve  on  these  roots, 
which  are  then  used  to  find  the  maxima  and  minima 
of  the  1st  derivative.  The  whole  process  is  repeated 
for  finding  the  roots  of  the  1st  derivative  and  then  for 
finding  the  roots  of  the  quartic  equation.  Redundant 
roots  are  thrown  out.  The  roots  are  then  arranged 
in  the  order 

Ds,  Df,  Us,  Uf 

i.e.,  downgoing  slow,  downgoing  fast,  upgoing  slow, 
upgoing  fast.  In  the  eventuality  that  the  slow 
mode  has  larger  attenuation  than  the  fast  mode 
then  the  slow  and  fast  modes  are  interchanged  in 
order  to  prevent  swamping. 
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EIGENV 


Computes  the  eigenvectors  associated  with  the  Booker 


quartic  roots .  For  any  given  quartic  root  the  charact¬ 
eristic  matrix  has  16  different  3x3  determinants 
which  may  be  obtained  by  elimination  of  the  ith 
row  and  jth  column.  Far  each  eigenvalue  the  two 
largest  of  the  16  determinants  are  found.  The 
eigenvectors  associated  with  each  of  the  two  largest 
determinants  are  then  calculated .  Thus  for  each 
eigenvalue  two  different  sets  of  equations  have  been 
used  to  calculate  two  eigenvectors .  The  eight  eigen¬ 
vectors  assoicated  with  the  four  eigenvalues  may  be 
combined  in  sixteen  differnet  ways  to  give  the  eigen¬ 
vector  matrix .  The  eigenvector  matrix  which  has  the 
maximum  determinant  is  chosen  because  the  columns 
of  that  particular  eigenvector  matrix  a  s  as  independent 
as  possible. 

BOBM  This  is  a  University  of  Pittsburgh  subroutine  which 

causes  the  computer  to  halt  at  the  next  M  pause" 
statement.  A  message  is  then  printed  to  the  operator 
to  mount  a  third  private  tape.  i.a. ,  tape  IW4. 

2.33  Program  for  Part  HT 

The  basic  function  of  Part  It!  is  to  perform  the  actual  Integration 
of  the  wave  equation.  The  forward  integration  proceeds  in  the  upward 
direction,  starting  deep  in  the  interior  of  the  Earth  and  proceeding 
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out  to  HMAX,  at  which  point  the  boundary  conditions  for  the 
disturbance  source  are  applied .  Because  of  the  problem  of  numerical 
swamping  (See  Paper  II  for  a  complete  discussion}  corrections 
frequently  have  to  be  made  to  the  solutions  in  the  various  horizontal 
slabs  in  order  to  preserve  good  linear  independence.  In  order  to  be 
efficient  the  corrections  are  all  made  in  one  backward  integration 
process,  which  is  performed  after  the  boundary  conditions  are  applied 
at  HMAX.  The  matrlcant  method  of  matrix  algebra  (Gantmacher,  1959; 
Pease,  1965;  Inoue  and  Horowitz,  1966b)  is  utilized  in  performing  the 
forward  integration. 

In  general  the  solution  to  the  wave  equation  is  given  by 

^ 

(taper  I,  equation  (14))  where  is  the  wave  field  solution  vector, 

0ff  the  4x4  integral  matrix  whose  columns  are  four  linearly  independent 
solutions  to  the  wave  equation,  and  the  superposition  vector'll consisting 
of  four  constants  given  by  the  boundary  conditions.  Because  the  dis¬ 
turbance  source  is  considered  to  be  above  the  layers  over  which  the 
integration  is  performed  two  of  the  superposition  constants  will  be  zero 
(See  Table  3 ,  taper  II) .  Consequently  only  two  independent  solutions 
are  necessary  for  the  final  solution.  These  will  correspond  to  the  incident 
slow  mode  and  the  incident  fast  mode  and  are  selected  by  the  variables  SBl 
and  SB2  in  the  computer  program  (See  Table  2.3S  below). 

The  following  basic  operations  performed  by  the  program  for 
Part  III; 
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1 .  The  synthesis  matrix  from  Part  II  is  fed  into  the  wave 
properties  matrix  (equation  (31),  Paper  II)  and  the  third  and  fourth 
columns  are  zeroed  out  because  only  two  independent  solutions  are 
required  as  a  result  of  the  boundary  conditions . 

2 .  The  average  propagation  tensor  is  found  by  linear  interpolation 
of  the  values  at  the  upper  and  lower  heights  of  a  slab. 

3 .  The  Booker  quartic  roots  and  the  hybridizing  matrices  (equations 
(19)  and  (27),  Paper  II)  are  calculated  in  accordance  with  whether  the 
slab  is  isotropic  or  anisotropic,  indicated  by  the  variable  IANISO 

which  was  given  a  value  in  Part  II. 

4.  The  modification  matrix  (equations  (9),  (10),  (18),  and  (22), 
Paper  II)  and  the  eigenvectors  are  then  calculated . 

5.  Undesirable  wave  components  are  excluded  from  the  wave 
properties  matrix  by  a  linear  combination  of  solutions  (section  5.2, 

Paper  II) . 

6 .  The  boundary  conditions  are  applied  at  HMAX  according  to 
whether  the  first  Independent  solution  (IRL  =  1)  or  the  second 
independent  solution  (IRL  =  2)  is  desired.  Either  case  utilizes  the 
same  forward  integration.  Only  the  backward  integration  is  different. 

7 .  On  the  backward  integration  the  correction  to  the  integral 
matrix  is  performed  by  use  of  the  CO'ft  matrix  (equation  (47),  Paper  II). 

The  amplitude  and  phase  factors  (equation  (32),  Paper  II)  and  the 
mixing  ratios  are  also  calculated.  The  mixing  ratios  are  defined  in  the 
following  manner. 
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e  =  >7*  t 

©  • 

ow 

where  *  is  the  principal  axis  transformation  of  the  propagation  tensor 
and  jW  contains  the  amplitude  characteristics .  Then 


where  f  is  the  column  matrix  of  elements  defined  as  the  mixing  ratios 
because  they  indicate  the  magnitude  of  the  "mixing"  of  the  characteristic 
waves,  which  are  the  columns  of  the  matrix. 

8 .  The  final  form  of  the  wave  fields  is  then  calculated  in  the 
background  integration. 

9 .  The  Booker  quartic  roots  <  the  electric  and  magnetic  wave 

field  magnitudes  and  the  mixing  ratios  for  the  first  and  second  independent 
solutions  are  stored  on  magnetic  tape  IW4  for  subsequent  automatic 
plotting . 

The  variables  read  into  the  computer  on  cards  as  data  are  as 
follows . 
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TABLE  2 . 35 


Variable 


HMAX 

HMIN 

v 

EUXB  (1) 
EUXB  (2) 
EUYB  (1) 
EUYB  (2) 


Description 

Upper  limit  (km.)  of  integration 
Lower  limit  (km.)  of  integration 
These  variables  represent  the  real  and  imaginary 
parts  of  the  x  and  y  components  of  the  electric 
field  of  downgoing  wave  at  the  upper  limit  HMAX 
of  integration,  i  .e . , 


tExL 


1 - 

J 

Kl 

where  all  quantities  are  evaluated  at  HMAX,  the 
subscripts  1  and  2  refer  to  the  1st  and  2nd  independent 
solutions,  and  the  superposition  constants  b’j  and  b'2 
are  given  by 


b’j  “  bj°jj 

where  the  Qjj  are  the  amplitude  and  phase  functions 
(equation  (32),  Paper  II)  and  the  bj  are  constants 
determined  by  the  boundary  conditions . 

There  are  two  different  ways  that  the  source  may 
be  specified  at  HMAX.  Either  the  constants  b'j  and 

wife 

b'2  are  given,  in  which  case  Eu  is  given  by  the 
equation  above,  or  else  Eu  is  specified  at  HMAX 
and  b'j  and  b^  are  then  predetermined  by  the  equations 
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EUXB  EUX2 


b2 

A 


EUYB  EUY2 

A 

EUX1  EUXB 

EUY1  EUYB 

A 


(EUX1)  (EUY2)  -  (EUY1)  (EUX2) 


When  it  is  preferred  that  b’j  and  be  specified 
then  EUXB  and  EUYB  are  given  the  value  zero  on  the 
input  data  card . 

SB1  (1)  These  variables  represent  the  real  and  imaginary 

SB1  (2) 

SB2  (1)  parts  of  b\  and  b'2,  respecti  Ay.  For  the  example 

SB2  (2) 

calculations  given  in  this  work ,  if  IRL  =  1  then 

b'_  =1.0  and  b'o  =0.0  and  if  IRL  =  2  then  b'  =0.0 
1  L  1 

and  b^  =  1.0.  Separate  backward  integrations  are 
performed  for  IRL  =  I  and  IRL  =  2. 

A  brief  description  of  the  various  subroutines  used  in  Part  III,  but 
not  in  Part  I  or  Part  II,  is  given  below. 

TABLE  2.36 

Subroutine  Description 

DECOMP  Used  to  check  the  linear  independency  oi  the  columns 

of  the  wave  properties  matrix.  Each  column  is 
normalized  with  respect  to  the  magnitude  of  that 
column.  It  is  then  easier  to  compare  different 
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columns  and  ascertain  the  degree  of  numerical 
swamping  and  degradation  in  the  linear  independence 
of  the  columns. 

LGP.ISO  Computes  the  hybridizing  matrices  for  the  isotropic 

case. 

LGRNG  Computes  the  hybridizing  matrices  for  the  anisotropic 

case. 

CMLOG  Calculates  the  natural  logarithm  of  a  complex  number. 

SKIP  This  is  a  University  of  Pittsburgh  subroutine  which 

takes  the  form  CALL  SKIP  (NF,  NR,  NTAPE) ,  where  NF 
is  the  number  of  file  marks  to  be  skipped,  NR  the 
number  of  records  to  be  skipped,  and  NTAPE  is  the 
number  of  the  tape  to  which  these  commands  apply. 


Copies  of  the  program  object  deck  and  flow  chart  may  be  obtained 
from  the  Air  Force  Cambridge  Research  Laboratories,  Ionospheric  Physics 
Labortory,  L.G.  Hanscom  Field,  Bedford,  Mass.  01730. 
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3.0  THE  TRANSMISSION  OF  ELF  WAVES  THROUGH  THE  IONOSPHERE 


3.1  Ion  Cyclotron  Resonance 

Figure  3.1  depicts  the  ion  cyclotron  resonance  profile  (i.e.,  the 
L  =  oc  surface  on  the  CMA  diagram)  for  the  particular  model  ionosphere 
used  in  this  work.  The  ions  primarily  responsible  for  each  portion  of  the 
curve  are  indicated.  That  part  of  the  profile  which  tails  off  to  the 
left  above  2000  km.  is  due  to  the  constancy  of  the  ion  mass  (H+)  and  the 

_3 

r  dependence  of  the  magnetic  field . 

In  this  analysis  a  two  component  plasma  was  assumed  and  the 
effects  of  ions  of  different  charge  to  mass  ratio  were  not  examined. 
Propagation  effects  in  a  plasma  containing  several  ion  species  have  been 
discussed  by  Stix  [1965].  Gurnett  et.  al.  [1965]  have  applied  the  theory 
of  propagation  in  a  multi-component  plasma  to  ion  cyclotron  whistlers 
and  Hayakawa  et.  al.  [1969]  have  calculated  the  dispersion  of  waves 
in  the  lower  exosphere  containing  multiple  ion  species  for  the 
frequency  range  0-150  c/s.  Also,  experimental  evidence  has  been 
found  for  the  existence  of  various  ion  resonances  in  the  ionosphere. 

Smith  et.  al.  [1964]  reported  the  discovery  of  the  proton  whistler  in 
satellites.  Barrington  et.  al.  [1966]  observed  the  helium  whistler  in 
Alouette  2  satellite  data.  More  recently  Gurnett  et.  al.  report  the 
observation  of  an  ion  cyclotron  whistler  due  to  an  ion  with  a  mass  to 
charge  ratio  of  8  amu/unit  charge.  This  ion  is  thought  to  be  either 
0++  or  Hejj,  but  its  identity  has  not  been  clearly  established.  There 
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is  every  reason  to  believe  that  other  ion  resonances ,  such  as 
No+,  O^,  N+,  0+,  etc.,  will  be  discovered . 

The  effect  that  multiple  ions  have  on  R  and  L  (the  squares  of 
the  indices  of  refraction  for  the  right  and  left  circularly  polarized 
waves,  respectively;  See  Stix  [1962])  at  1000  km  for  the  model 
ionosphere  used  iu  this  work  may  be  seen  in  Figure  3.2.  In 
this  graph  collisions  are  excluded  and  L  goes  to  infinty  at  the 
gyrofrequency  for  each  ion  present  in  the  plasma.  For  a  two  component 
plasma  there  will  appear  only  one  resonance,  which  corresponds  to 
the  gyrofrequency  of  the  effective  ion.  In  addition  between  each 
gyrofrequency  there  appears  a  point  at  which  D=  1/2  (R  -  L)  =  o,  ie,  R  and 
L  cross  over  each  other.  The  frequency  at  which  this  occurs  is  called  the 
"crossover"  frequency.  Smith  and  Brice  [1964]  find  that  both  modes  are 
linearly  polarized  at  the  crossover  frequency.  Polarization  reversal 
occurs  for  a  wave  passing  through  the  crossover  frequency.  This 
situation  corresponds  to  the  phase  velocity  surfaces  being  tangent  in 
the  direction  parallel  to  the  magnetic  field .  At  this  point  the  R  and  L 
labels  are  interchanged .  Such  a  phenomena  does  not  occur  for  a  two 
component  plasma  because  D  can  strictly  vanish  only  for  a  plasma 
containing  positive  ions  of  different  charge-to-mass  ratio,  provided  the 
magnetic  field  and  electron  density  are  not  zero  (See  problem  2, 
chapter  1 ,  Stix  [1962]) .  In  fact  D  may  vanish  only  at  frequencies 
greater  than  the  m  inimum  ion  cyclotron  frequency  and  less  than  the 
maximum  ion  cyclotron  frequency. 
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The  discussion  so  far  has  been  limited  to  the  case  of  a  collisionless 


plasma.  Including  collis'^ns  causes  the  quantities  R,  L,  and  D  to  contain 
imaginary  as  well  as  real  parts .  A  plot  of  the  real  parts  of  R  and  L  with 
and  without  collisions  for  a  model  ionosphere  with  multiple  ions  and  also 
for  the  case  of  an  effective  ion  suffering  collisions  is  shown  in 
Figure  3.3.  A  representative  altitude  of  200  km  is  chosen .  The 
subscripts  "O"  and  "i"  refer  to  the  quantities  calculated  excluding 
collisions  and  including  collisions,  respectively,  but  for  a  multi- 
component  plasma  in  both  cases.  The  subscript  "effM  refers  to  a  two 
component  plasma  containing  electrons  and  an  effective  ion,  with 
both  species  suffering  collisions. 

It  is  evident  from  Figure  3 . 3  that  the  Inclusion  of  collisions  at 
200  km  has  no  significant  affect  on  the  value  of  R.  This  is  consistent 
with  the  result  obtained  by  Gumett  et.  al.  [1965]  that  the  real  part 
of  R  for  the  collisionless  case  is  the  same  as  the  real  part  of  R  with 
collisions  when  R  is  expanded  to  first  order  in  the  parameter  •  Oe . , 
the  ratio  of  collision  frequency  to  angular  wave  frequency) .  This  same 
result  is  seen  to  be  true  for  L  when  sufficiently  far  away  from  the 
gyrofrequencles .  However,  collisions  do  significantly  affect  the 
/slue  of  L  near  the  gyrofrequencles .  The  value  of  L  is  moderated  by 
collisions  and  no  longer  tends  to  infinity  at  the  gyrofrequencles .  Alto  L 
may  go  through  zero  near  the  gyrofrequencles  and  at  some  point  between 
any  two  ion  gyrofrequencles  may  become  equal  to  R.  Thus  even  with 
collisions  included  there  may  still  be  frequencies  between  the  ion 


gyrofrequencies  at  whic’  the  real  part  of  D  is  zero. 

The  question  still  remains  as  to  whether  the  vanisning  of  the  r^al 
part  of  D  corresponds  to  a  polarization  reversal.  Gurnet  et.  al.  claim 
that  when  the  magnitudes  of  the  imaginary  parts  of  R,L,  P,  and  H 
(H  =  RL  -  PS)  are  much  less  than  the  real  parts  of  the  corresponding 
quantities  then  there  exists  a  critical  coupling  angle  Oc  near  DQ=  o 
determined  by  the  relation 

where  the  subscripts  ''o"  and"!"  refer  to  the  real  and  imaginary  parts, 
respectively.  Thus,  in  the  presence  of  collisions  the  wave  normal 
surfaces  at  »  o  intersect  over  a  cone  of  directions  making  an  angle 
9C  with  the  magnetic  field.  This  cone  is  called  the  critical  coupling 
cone.  If  a  wave  passes  through  this  critical  coupling  cone  near  DQ  *  o 
a  wave  of  opposite  polarization  is  generated . 

However,  in  Stix'  notation,  the  condition  of  critical  coupling  is  the 
2 

vanishing  of  F  ,  i .  e . , 

F  *  »o 

Then,  in  the  case  of  small  collisions,  the  exact  condition  for  critical 
coupling  is 
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and  the  critical  coupling  angle  ^  is  given  by 


,  ,  i  Iw  1 

<•  =  *  ^K*ou.-psi  J 


=y  ,  *eCRL-Psy 

Therefore  critical  coupling  does  not  necessarily  take  place  near  the  level  of 


occur. 


RttP>o.  In  fact  even  if  \>  critical  coupling  may 

Also,  with  collisions  included,  polarization  reversal  may  arise  from  the 

2 

quadrant  change  of  F  in  the  complex  plane  at  a  certain  angle  without 
causing  critical  coupling. 

Consideration  of  an  effective  ion  and  collisions  gives  rise  to  a 
moderated  resonance  effect  only  in  the  vicinity  of  the  gyrofrequency 
corresponding  to  the  effective  ion  >ee  Figure  3.3).  Also,  there  may 
exist  one  crossover  frequency  (ie. ,  DQ  =  o)  when  only  the  effective 
ion  is  considered,  provided  collisions  are  included.  However,  at  most 
one  such  frequency  exists. 

Therefore,  this  discussion  serves  to  show  that  even  in  the  presence 
of  collisons  it  is  necessary  to  include  all  the  ions  in  order  to  determine 
the  wave  characteristics  of  a  disturbance  propagating  with  a  wave  frequency 
near  any  of  the  ion  gyrofrequencies . 


1-33 


3.2  Localized  disturbance  at  100  c/s 

3.21  Booker  Quartic  Roots 

T^ble  3.1  summarizes  the  nature  of  the  slow  and  fast  waves  at 
100  c/s  in  the  various  regions  of  the  Earth's  environment.  Figure  3.4 
shows  the  altitude  dependence  of  the  Booker  quartic  roots  for  the 
slow(s)  and  fast(F)  waves.  The  abbreviations  QR  and  QI  stand  for  the 
real  and  imaginary  parts  of  the  quartic  roots,  respectively.  The 
indices  of  refraction  for  propagation  in  the  z  -  direction  are  obtained 
by  multiplying  the  Q  values  by  a  factor  of  .075.  Downgoing  and 
uogoing  waves  are  represented  by  the  same  curves  and  differ  only  by 
a  minus  sign. 

The  Booker  quartic  roots  are  almost  entirely  real  in  the  region 
above  1070  km.  The  slow  and  fast  waves  are  quite  similar  in  this 
region  and  are  both  propagating.  At  1070  km  the  ion  cyclotron 
resonance  is  encountered.  Below  this  level  the  slow  and  fast  waves 
are  interchanged  and  the  fast  wave  is  evanescent.  This  resonance 
condition  corresponds  to  a  crossing  of  the  L  =  °o  surface  in  parameter 
space  on  the  CMA  diagram.  The  magneto-ionic  quantity  Yj  is  equal 
to  one .  As  we  shall  see  later  from  the  mixing  ratios  and  coupling 
coefficients  this  resonance  causes  the  L  mode  to  be  almost  completely 
ieflected  back  to  the  magnetosphere.  The  slow  mode  is  highly 
propagating  down  to  the  neutral  atmosphere.  Above  the 
resonance  level  the  ionosphere  is  slightly  anisotropic.  Below  the 
resonance  level  the  ionosphere  is  highly  anisotropic. 
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PROPAGATION  CHARACTERISTICS 


TABLE  3.1  Medium  characteristics  at  100  c/s  for  a  localized 
disturbance  incident  at  an  oblique  angle  from  the  vertical  at  5000  Km. 
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Figure  3.4 


In  the  neutral  atmosphere  the  quartic  roots  are  pure  imaginary 
and  both  modes  are  evanescent.  This  effect  is  due  to  the  presence  of 
a  horizontal  component  in  the  propagation  vector.  In  the  neutral  atmosphere 
the  electron  and  ion  densities  are  zero.  Therefore,  in  Stix'  notation: 

r  =  L=  P=  S=1 
A  =  C  =  1 
B  =  2 

The  dispersion  relation  then  becomes 

d  1 

/I  -  in  +i  =  o 

or  Cn-0  =  O 

*  =  I  =  + 

thus  <  =  i-  fa 

t  x  *  * 

If  oT«c,  then 

x  «x 

Y\  <  O 
£ 

and  hence  /I  is  pure  imaginary. 

Thus  the  neutral  atmosphere  acts  as  an  elastic  medium  for  the  vertically 
propagating  component  of  the  particular  localized  disturbance  being  con¬ 
sidered  .  If  /\£  had  been  greater  than  c  then  o  <  |  and  the 

neutral  atmosphere  would  have  been  dispersive  for  the  vertically 
propagating  component. 

Inside  the  Earth  the  magnitudes  of  the  real  and  imaginary  parts 
of  the  quartic  roots  are  all  equal.  Thus  the  Earth's  interior  is 
isotropic  and  "dense  resistive",  ie.  behaves  like  a  metal. 

3.22  Magneto-ionic  Coupling  Coefficients 

Figure  3.5  shows  the  magneto-ionic  coupling  coefficients  (See 
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MAGNETO  -  IONIC  COUPLING 


Inoue  and  Horowitz,  1966a)  for  the  localized  disturbance  at  100  c/s. 

The  abbreviations  S,  F,  M,  R  stand  for  slow,  fast,  mode  coupling  and 
reflection,  respectively.  Thus  RS  means  reflection  of  the  slow  wave. 

M  means  generation  of  a  downgoing  slow  wave  from  a  downgoing  fast  wave 
or  an  upgoing  fast  wave  from  an  upgoing  slow  wave,  etc.  MR  means 
generation  of  a  upgoing  fast  wave  from  a  downgoing  slow  wave  or  an 
upgoing  slow  wave  from  a  downgoing  fast  wave,  etc. 

Reflection  is  the  dominant  type  of  coupling  at  almost  all  heights . 

At  1070  km  reflection  type  coupling  for  the  slow  wave  becomes  very  large 
because  of  the  ion  cyclotron  resonance.  For  the  fast  wave  below  this 
level  the  magnitude  of  the  coupling  is  larger  than  or  comparable  to  the 
quartic  roots .  For  the  slow  wave  below  the  resonance  level  coupling  is 
of  less  importance  than  the  quartic  rods  down  to  the  D  region.  In 
the  lower  D  region  coupling  becomes  very  large  for  the  slow  and  fast 
modes .  Mode  reflection  type  coupling  also  appears .  The  plane  wave 
case  produces  no  mode  reflection  type  coupling  in  the  lower  D  region. 
Therefore  this  effect  is  due  to  the  localized  nature  of  the  disturbance. 

Reflection  type  coupling  occurs  within  the  surface  layer  of  the  earth 
and  is  a  maximum  at  a  depth  of  about  2.5  km.  This  reflection  region 
corresponds  to  the  rapid  change  in  the  earth  conductivity.  In  this  region 
the  coupling  coefficients  are  in  general  smaller  than  the  quartic  roots 
but  are  of  comparable  magnitude  at  a  depth  of  2.5  km. 
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3.23  Characteristic  Wave  Mixing  Ratios 

The  mixing  ratios  (See  section  2.33  for  definition)  for  the 
incident  fast  wave  are  given  in  Figure  3.6.  Above  1070  km.  the  resultant 
wave  fields  consist  primarily  of  a  downgoing  fast  (DF)  wave  but  with 
some  partial  reflection  occurring  as  may  be  seen  by  the  presence  of  a 
small  amount  of  the  upgoing  fast  wave  (UF) .  The  mixing  ratio  for  the 
upgoing  fast  wave  exhibits  a  standing  wave  nature  due  to  reflection  from 
the  ion  cyclotron  resonance  level.  The  downgoing  slow  wave  (DS) 
is  generated  by  mode  coupling,  which  may  also  be  seen  in  the  coupling 
coefficients  in  Figure  3.5.  The  upgoing  slow  (US)  wave  is  generated 
by  mode  reflection  type  coupling,  again  also  seen  in  Figure  3.5. 

However,  these  last  two  waves  offer  only  a  small  contribution  to  the 
total  wave  field . 

From  1070  km.  down  through  the  D  region  the  dominant  mixing 
ratio  is  that  of  the  riowngoing  slow  wave .  (Remember  that  below  the 
critical  reflection  layer  the  R  mode  becomes  a  slow  wave) .  Some 
partial  reflection  is  seen  through'''”':  this  altitude  range  as  evidenced 
by  the  mixing  ratio  for  the  upgoing  slow  wave.  Between  700  and 
1000  km.  there  is  a  standing  wave  due  to  the  ion  mass  changing  very 
rapidly.  The  standing  wave  between  300  and  600  km.  is  caused  by  the 
upper  slope  of  the  F  maximum  electron  density.  There  is  also  a 
standing  wave  centered  on  the  F  region  maximum. 

In  the  D  region  the  mixing  ratios  for  the  downgoing  and  upgoing 
fast  waves  grow  rapidly  as  a  consequence  of  the  mode  reflection  type 
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coupling  already  seen  in  Figure  3.5.  This  effect  is  due  to  the 
localixed  nature  of  the  disturbance . 

In  the  neutral  atmosphere  and  in  the  ground  the  transverse  magnetic 
(TM)  wave  is  larger  than  the  transverse  electric  (TE).  The  earth  acts  as 
a  good  conductor  and  almost  complete  reflection  occurs  as  seen  from 
the  approximate  equality  of  the  downgoing  and  upgoing  parts  of  the 
TM  and  of  the  TE  waves.  The  mixing  ratios  are  attenuated  logarithmically 
in  the  ground  and  only  downgoing  waves  are  present.  Although  the 
coupling  coefficients  indicate  that  reflection  occurs  for  both  the  transverse 
magnetic  and  transverse  electric  modes  between  about  1  and  3  km .  depth 
this  process  has  negligible  effect  on  the  resultant  mixing  ratios  and 
wave  fields  due  to  the  fact  that  the  quartic  roots  are  larger  than  the  coupling 
coefficients  and  are  therefore  more  important  in  determining  the  final 
wave  field  characteristics . 

The  mixing  ratios  for  the  incident  slow  wave  may  be  seen  in 
Figure  3.7.  The  equality  of  the  mixing  ratios  for  the  downgoing  slow  and 
upgoing  slow  waves  indicates  that  complete  reflection  takes  place  at 
1070  km, ,  the  ion  cyclotron  resonance  level.  Thus  the  slow  mode 
(ie.  the  L  mode)  is  completely  non- penetrating  below  this  level.  This 
result  is  in  agreement  with  the  prediction  of  the  CMA  diagram  that 
the  wave  normal  surface  for  the  L  mode  disappears  upon  crossing  the 
L-oo  surface  from  lower  to  higher  frequencies. 
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3.24  Wave  Fields 


The  wave  fields  of  the  incident  fast  wave  are  shown  in  Figure  3.8. 

The  magnitudes  of  the  electric  and  magnetic  fields  are  ail  relative  to 
the  values  at  4000  km.  If  the  magnetic  field  values  are  multiplied  by 
a  factor  of  0.25  the  B  field  is  then  obtained  in  gammas  (1  gamma  = 

10  °  gauss).  The  electric  field  is  read  directly  in  volts/meter.  The 
letters  x,  y,  z  stand  for  the  corresponding  components  of  the  field 
quantities . 

The  wave  fields  exhibit  the  same  standing  wave  characteristics 
as  the  mixing  ratios  in  the  F  region.  The  standing  wave  above  1070  km. 
is  due  to  partial  reflection  above  the  ion  cyclotron  resonance  level.  The 
standing  wave  between  300  and  600  km.  is  caused  by  partial  reflection 
from  the  upper  slope  of  the  F  maximum  electron  density.  A  startdinq  wave 
Is  also  centered  on  the  F  layer  maximum.  Standing  waves  are  visible 
in  the  wave  fields  dvwn  through  the  E  region  but  are  less  obvious  in  the 
D  region  and  neutral  atmosphere.  In  the  D  region  mode  coupling  produces 
the  separation  in  the  x  and  y  components  of  the  electric  and  magnetic 
fields  and  an  increase  in  the  z  -  component  of  the  electric  field.  The 
e' metric  field  remains  essentially  vertical  throughout  the  neutral  atmosphere 
and  must  cause  charge  separation  to  occur  between  the  Earth's  surface 
and  the  lower  ionosphere.  In  the  ground  the  fields  are  attenuated 
logarithmically. 

The  wave  fields  for  the  incident  slow  wave  display  a  standing 
wave  effect  above  1070  km.  due  to  reflection  from  the  ion  cyclotron 
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resonance  level .  The  fields  are  completely  zero  below  this  level  and 
further  substantiate  the  fact  that  the  L  mode  is  non- pentra ting  beyond 
the  ion  resonam  i  level. 

3.25  Wave  Field  Polarization 

The  polarizations  of  the  electric  and  magnetic  wave  fields  for 
the  incident  fast  wave  are  given  in  Figure  3.9.  The  labeling  at 
the  top  of  each  diagram  specifies  the  plane  in  which  the  polarization  is 
plotted.  For  example,  "YZ"  implies  that  the  y  component  of  the  field 
quantity  is  plotted  along  the  horizontal  direction  and  the  z  component 
is  plotted  along  the  vertical  direction.  The  number  under  each  polarization 
curve  represents  the  altitude  in  km.  for  which  that  polarization  plot 
is  valid.  The  Earth's  magnetic  field  is  directed  inward  and  rotation 
in  the  clockwise  direction  represents  right  handed  polarization.  To 
obtain  the  direction  of  rotation  the  gap  in  the  polarization  curve  must 
first  be  found.  That  side  of  the  gap  which  Is  marked  with  an  "x" 
represents  the  beginning  point.  At  90°  from  the  "x"  is  a  dot  and  the 
rotation  sense  is  obtained  by  proceeding  from  the  "x"  to  the  dot. 

At  each  altitude  the  plot  is  normalized  to  the  larger  of  the  two 
components  being  plotted .  Thus  in  order  to  determine  the  relative 
sizes  of  plots  at  different  altitudes  one  must  determine  the  wave 
field  magnitudes  from  Figure  3.8. 

In  the  x-y  plane  the  electric  and  magnetic  fields  are  right 
circularly  polarized  down  to  the  D  region.  Due  to  the  mode  coupling 
which  produces  a  separation  in  the  x  and  y  components  of  the  electric 
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3.9 


and  magnetic  fields  in  the  D  region  the  fields  become  elliptically 
polarized  in  the  lower  ionosphere.  At  50  km.  the  electric  field  becomes 
left  elliptically  polarized  and  then  changes  back  to  right  elliptically 
polarized  at  ground  level .  Thus  the  fast  mode  which  is  polarized  in 
the  right  hand  sense  at  4000  km.  is  also  right  handed  polarized  at 
ground  surface. 

However,  on  the  Earth's  surface  the  x  and  y  components  of 
the  electric  field  are  smaller  than  the  z  component  by  three  and 
four  orders  of  magnitude  respectively.  Thus  the  xz  and  yz  planes 
ate  more  important  than  the  xy  plane.  In  these  planes  the  electric  field 
on  the  ground  is  linearly  polarized. 

The  incident  slow  mode  is  left  circularly  polarized  throughout  the 
region  above  the  ion  cyclotron  resonance  level. 

5.3  Comparison  of  Plane  Wave  with  Localized  Disturbance  at 
i  00  c/s 

Comparison  of  the  effects  of  the  plane  wave  and  localized  dis¬ 
turbance  shows  very  little  difference  between  the  incident  slow  modes, 
doth  are  non -penetrating  below  about  1070  km.  and  in  both  cases  the 
downgoing  slow  wave  is  reflected  in  the  form  of  an  upgoing  slow  wave. 
The  major  difference  seems  to  be  that  the  wave  fields  of  the  plane 
wove  are  TT  out  of  phase  with  those  of  the  localized  disturbance. 
However,  comparison  of  the  incident  fast  modes  does  reveal  some 
differences.  Because  there  is  no  propagation  in  the  horizontal  direction 
the  plane  wave  propagates  vertically  through  the  neutral  atmosphere 
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with  a  phase  velocity  equal  to  the  speed  of  light.  Thus  the  neutral 
atmosphere  behaves  like  free  space  for  the  plane  wave  but  appears  as 
an  elastic  evanescent  region  to  the  localized  disturbance. 

Also,  mode  coupling  and  mode  reflection  type  coupling  appear 
for  the  localized  disturbance  in  the  lower  D  region  but  no  such  coupling 
appears  for  the  plane  wave.  Reflection  of  the  slow  mode  however 
still  occurs  throughout  the  ionosphere  for  the  plane  wave  just  as  it  does 
when  the  localized  disturbance  effect  is  introduced.  The  z  components 
of  the  electric  and  magnetic  fields  are  zero  everywhere  for  the  plane 
wave  and  the  electric  field  polarization  in  the  xy  plane  remains  right 
circularly  polarized  throughout  the  ionosphere  and  the  neutral 
atmosphere.  In  contrast  the  localized  effect  produces  a  z  component  of 
the  magnetic  field  throughout  the  ionosphere  and  atmosphere  and  a  z 
component  of  the  electric  field  in  the  lower  D  region  and  neutral 
atmosphere.  The  z  components  are  due  to  the  additive  interference  of 
the  upgoing  and  downgoing  waves .  Therefore  only  the  localized 
disturbance  gives  rise  to  charge  separation  between  ground  level  and 
the  lower  ionosphere.  Also  the  localized  disturbance  wave  field 
polarizations  changed  considerably  in  the  lower  D  region  due  to 
strong  coupling. 

Within  the  Earth  the  x  and  y  components  for  both  the  electric 
and  magnetic  fields  are  equal  in  the  plane  wave  case  but  the  localized 
disturbance  gives  rise  to  an  inequality  in  those  components. 
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3.4  Ionospheric  Electric  and  Magnetic  Field  Measurements 

Very  few  field  measurements  exist  for  the  ELF  range  around 
100  c/s  but  several  experiments  have  been  performed  [Maynard  et.  al. , 
1369;  Shawhan  et  al.,  1968;  Scarf  et.  al.,  1968;  Gurnett  et.  al . , 

1969].  Various  difficulties  have  been  encountered  and  not  all 
measurements  are  in  agreement.  Maynard  et.  al.  measured  the 
complete  wave  characteristics  of  several  whistlers,  whose  source 
was  the  local  ionosphere  waveguide,  between  approximately  100  and 
200  km.  with  a  ;e-Tomahawk  sounding  rocket.  In  particular,  one 
whistler  measured  at  171.2  km.  was  elliptically  polarized  and  had 
the  lowest  frequency  component  observed,  ie.  Go  c/s,  with  an 
electric  field  value  between  100  and  200  uv/'tn.  Shawhan  et.  al . 
measured  electric:  and  magnetic  fields  with  the  Javelin  8.45  sounding 
rockot.  However,  the  measurements  were  rather  broadband 
(30-650  c/s)  and  indicated  electrostatic  noise  (ie.  no  magnetic  field;. 
Also  Ehawhan  et.  ai.  were  not  sure  if  these  measurements  were 
indicative  of  a  natural  phenomena  or  due  to  payload  interaction  with  the 
ionosphere.  Scart  et.  al .  performed  electric  and  magnetic  field 
measurements  with,  sensing  equipment  on  board  the  OV3-3  satellite. 
Most  observations  were  electrostatic  but  one  case  was  electromagnetic. 
At  an  altitude  of  364  km.  the  maximum  magnetic  field  strength  was 
80  milligammas  at  a  frequency  of  525  c/s.  The  electric  field 
strength  was  on  the  order  of  100  microvolts/meter .  Low  frequency 
measurements  were  tampered  by  spacecraft  noise.  Gurnett  et.  al . 
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obtained  observations  of  electric  and  magnetic  fields  with  the  Injun 
5  satellite,  launched  on  August  8,  1968,  into  an  elliptical  polar  orbit 
with  an  apogee  of  2528  km.  and  a  perigee  of  677  km.  The  Injun  5 
experiment  had  the  capability  of  measuring  electric  and  magnetic  fields 
of  plasma  wave  phenomena  over  two  frequency  bands,  0.3  to 
10  kc/s  and  30  to  650  c/s.  The  maximum  magnetic  field  strenths  for 
the  latter  band  were  about  50  miliigammas.  The  maximum  electric 
field  strengths  were  on  the  order  of  3.5  mv/m  for  the  lower 
frequency  band  but  in  general  the  wide-band  electric  field  strengths 
are  much  less,  of  the  order  of  100  uv/m  or  less,  ~hese  electric  field 
strengths  are  smaller  than  those  reported  by  Scarf  et.  al .  Gurnett 
et.  al.  found  that  the  electric  and  magnetic  field  strengths  vary  together, 
suggesting  that  the  field  variations  are  due  to  electromagnetic  waves, 
possibly  ELF  hiss  where  the  variations  are  smooth,  and  discrete  emissions 
(chorus)  where  the  fluctuations  are  rapid. 

A  typical  value  of  100  miliigammas  is  assumed  for  the  magnetic 
wave  field  strength  at  360  km.  and  is  used  in  conjunction  with  the  full 
wave  results  to  predict  electric  anu  magnetic  wave  field  strengths  at 
other  altitudes  (See  Table  5.7,  Ghapt-r  :>). 
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4.0  THE  TRANSMISSION  OF  UI.F  WAVES  THROUGH  THE  IONOSPHERE 


4.1  Localized  Disturbance  at  3.0  c/s 

4.11  Cooker  Quartic  Roots 

Table  4.1  summarizes  the  properties  of  the  slow  and  fast  wave 
disturbances  at  3.0  c/s  in  the  ionosphere,  neutral  atmosphere  and  the 
Earth’s  interior.  The  Booker  quartic  roots  are  presented  in  Figure  4.1. 
The  labeling  used  in  this  figure  is  the  same  as  that  explained  >n  section 
3.21.  The  actual  values  of  the  index  of  refraction  for  propagation 
in  the  vertical  direction  are  obtained  by  multiplying  the  values  read  from 
the  curves  in  Figure  4 . 1  by  a  factor  of  2  .SO  . 

Comparison  of  the  real  and  imaginaiy  parts  of  the  quartic  roots 
reveals  that  both  the  slow  and  fast  modes  are  propagating  down  to 
about  130  km.  The  region  above  1000  km.  is  quosi-isotrop; ; .  ,V 
1000  km.  mode  switching  occurs  and  the  slow  .vave  t  , -comes  a 
fa  '*  wave  a’  d  the  fust  w.vo  becomes  a  slow  wave .  At  1-0 
km.  mode  switching  ■*'  din  occurs  with  the*  slow  mode  changing  back 
to  the  fast  mode  and  vice  versa.  In  section  4  .  15  wo  shall  see 


that  the  wave  Helds  undergo  a  polarization  reversal  at  1 00 u  km.  I  ut 
not  at  l  50  k;  > .  This  region  between  lUftO  km.  and  1  >0  km.  is  essentially 
isotropic,  between  220  km.  and  U»0  km.  the  rea*  *  arts  <v  the  s!  >w 
and  fast  wn  .*e  g  urt.se  roots  are  mthsur  : .  ishai  le  } lowe-.-er .  the  fa  r' 

wove  is  evanescent  between  IS.  km-  and  the  Earth's  surface.  The  slew 
wave  is  propagating  down  io  about  70  km.  at  which  it  becomes 
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'.aracteristics  at  3.0  c/s  for  a  localized 
!  oblique  angle  from  the  vertical  at  5000  Km. 


evanescent  and  remains  as  such  down  to  ground  level.  The  lower  D 
region  between  about  70  km.  and  50  km.  is  quasi-isotropic.  The 
phase  velocities  of  the  fast  and  slow  waves  reach  their  minimum  values 
around  250  km,  and  120  km.,  respectively . 

The  neutral  atmosphere  is  isotropic  and  the  quartlc  roots  for 
this  region  are  pure  imaginary  due  to  the  localization  of  the 
disturbance  (See  section  3.21).  The  interior  of  the  Earth  is  also  isotropic 
and  behaves  as  a  metal. 

4.131  Effect  of  Charge  Exchange  on  the  Booker  Quartlc  Roots 
Comparison  of  the  quartic  roots  obtained  by  including  charge 
exchange  with  the  quartic  roots  obtained  by  excluding  charge  exchange 
collisions  is  shown  in  Figure  4.2.  The  real  parts  of  the  quartic  roots 
are  essentially  unaffected  by  charge  exchange  between  150  and  800  km. 
altitude.  On  the  other  hand  the  imaginary  parts  of  the  Booker  quartic 
roots  are  enhanced  by  as  much  as  an  order  of  magnitude.  In  the  altitude 
range  150  to  250  km.  the  imaginary  parts  of  the  quartic  roots  with  charge 
exchange  included  are  the  same  order  of  magnitude  as  the  real  parts. 
Hence  inclusion  of  charge  exchange  collisions  significantly  enhances 
*he  absorption  of  the  disturbance  in  the  region  3  50-250  km.,  where 
charge  xchange  ofO+  with  O  is  the  dominant  process.  Figure  4.5  and 
4.6  show  that  the  electric  and  magnetic  wave  fields  are  indeed  greatly 
reduced  in  amplitude  in  this  region. 

4.12  Magneto-ionic  Coupling  Coefficients 

Figure  4 . 3  reveals  the  altitude  dependence  of  the  magneto- 
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MAGNETO  -  IONIC  COUPLING 

Figure  4 . 3 


ionic  coupling  coefficients  for  the  localize!,  isturbance  at  3.0  c/s. 

These  curves  have  the  same  nomenclature  attached  to  them  as  that 
explained  in  section  3.22. 

Reflection  type  coupling  is  dominant  down  to  about  1000  km. 

At  this  altitude  there  is  a  strong  mode  coupling  layer  caused  by  the 
equality  of  the  fast  and  slow  wave  quartic  roots,  ie .  the  medium  rapidly 
approaches  isotropic  conditions.  Reflection  type  coupling  for  both 
the  slow  and  fast  modes  once  again  becomes  the  major  coupling  process 
below  about  900  km.  and  remains  dominant  down  to  the  lower  D  region. 

In  the  F  region  mode  coupling  and  mode  reflection  type  coupling  are 
also  present  due  to  the  quasi-isotropic  nature  of  the  F  region.  These 
coupling  processes  become  small  in  the  E  region  because  the  medium 
is  anisotropic  and  reappear  in  the  lower  D  region  due  to  the  localized 
disturbance  (See  section  3.22).  Thus  coupling  is  of  importance  in 
the  D  and  F  regions  but  not  in  the  E  region.  The  quartic  roots  are  larger 
than  the  coupling  coefficients  throughout  the  F  region  except  for  the  region 
around  1000  km.  where  mode  coupling  becomes  very  large.  The  fast 
wave  quartic  roots  are  smaller  than  the  coupling  coefficients 
throughout  the  D  and  E  regions .  The  coupling  coefficients  become 
larger  than  the  slow  wave  quartic  roots  below  about  95  km.  Coupling 
in  the  neutral  atmosphere  is,  of  course,  zero. 

As  with  the  100  c/s  case,  reflection  type  coupling  occurs  in 
the  surface  layer  of  the  Earth  and  reaches  a  maximum  at  a  depth  of 
about  2.5  km.  This  reflection  process  corresponds  to  the  rapjd 
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change  in  the  ground  conductivity.  The  coupling  coefficients  in  this 
region  are  smaller  than  the  quartic  roots  but  become  of  similar  magnitude 
around  a  depth  of  2.5  km.  The  medium  inside  the  Earth  is  isotropic 
and  behaves  as  a  metal. 

4.13  Characteristic  Wave  Mixing  Ratios 

The  characteristic  wave  mixing  ratios  for  the  incident  slow 
and  fast  waves  at  3.0  c/s  are  depicted  in  Figure  4.4.  The  labeling 
used  for  these  curves  is  explained  in  section  3.23. 

In  the  region  above  1000  km.  the  resultant  wave  for  the  incident 
slow  mode  consists  primarily  of  the  downgoing  slow  wave.  The  presence 
of  the  upgoing  slow  wave  is  indicative  of  reflection  taking  place 
throughout  the  region.  The  downgoing  and  upgoing  fast  waves  are  also 
present  and  are  generated  by  mode  coupling  and  mode  reflection  type 
coupling.  In  the  F  region  below  1000  km.  the  downgoing  fast  and 
downgoing  slow  waves  are  nearly  equal  in  magnitude  due  to  the  strong 
coupling  associated  with  this  quasi-isotropic  region.  The  upgoing  fast 
and  upgoing  slow  waves  are  also  present  due  to  reflection  from  the  F 
region  maximum.  The  standing  wave  formation  in  these  upgoing  waves 
is  quite  obvious.  Around  200  km.  the  real  and  imaginary  parts  of  the 
quartic  roots  are  equal  as  in  a  conductor.  Therefore  this  level  will 
behave  as  the  upper  conducting  wall  of  a  cavity  with  the  other  conducting 
wall  being  the  Earth's  surface.  This  fact  is  borne  out  when  we  look  at 
the  standing  wave  formation  observed  in  the  wave  fields  in  the  next 
section.  Below  this  level  in  the  D  and  E  regions  the  downgoing  slow  and 
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upgoing  slow  waves  are  almost  equal,  indicating  that  reflection  is 
continuously  taking  place.  In  the  neutral  atmosphere  the  incident 
slow  wave  is  comprised  of  both  transverse  electric  and  transverse 
magnetic  waves  with  the  transverse  electric  beiny  dominant.  Both  the 
TE  and  TM  modes  are  reflected  quite  well  at  ground  surface. 

Reflection  of  the  transverse  electric  wave  is  also  een  to  take  place 
within  the  Earth  and  the  maximum  of  the  reflected  wave  occurs  at  a 
depth  of  about  two  km.,  near  the  conductivity  maximum. 

For  the  incident  fast  wave  above  1000  km.  all  of  the 
characteristic  waves  are  of  the  same  order  of  magnitude.  The  upgoing 
waves  are  generated  by  reflection  and  mode  reflection  type  coupling. 
Although  this  region  is  anisotropic  it  is  sufficiently  close  to  being 
isotropic  that  mode  coupling  and  mode  reflection  type  coupling  are 
important.  In  fact  the  downgoing  slow  wave  is  at  times  larger  than  the 
downgoing  fast  wave.  In  the  F  region  below  1000  km. ,  where  the  slow 
wave  becomes  a  fast  wave  and  vice  versa,  the  downgoing  slow  wave  is 
dominant  and  a  standing  w  /e  formation  occurs  in  the  reflected  upgoing 
slow  wave.  The  downgoing  fast  and  upgoing  fast  waves  are  also  present 
due  to  mode  coupling  and  mode  reflection  type  coupling,  respectively, 
resulting  from  the  quasi-isotropic  nature  of  this  region.  In  the  D  and  E 
regions  the  reflected  upgoing  fast  wave  is  smaller  than  the  downgoing 
and  upgoing  slow  waves,  which  are  generated  by  mode  coupling  and  mode 
reflection  type  coupling,  respectively.  Thus  the  D  and  E  regions  will 
not  act  to  store  energy  for  the  incident  fast  wave  as  readily  as  for  the 
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incident  slow  wave.  This  fact  will  be  easily  seen  from  inspection  of 
the  wave  fields.  In  the  neutral  atmosphere  the  incident  fast  wave  is 
transverse  electric  and  is  reflected  quite  well  at  the  Earth's  surface. 

Inside  the  Earth  the  incident  fast  wave  is  primarily  comprised  of  the 
downgoing  transverse  electric  mode.  However  reflection  produces  a 
small  amount  of  the  upgoing  transverse  electric  mode  which  reaches 
a  maximum  at  a  depth  of  about  two  km . ,  near  the  conductivity 
maximum . 

4.14  Wave  Fields 

The  wave  fields  of  the  incident  slow  wave  are  shown  in 
Figure  4.5.  The  magnitudes  of  the  electric  and  magnetic  fields 
are  all  relative  to  the  values  at  4000  km.  If  the  magnetic  field  values 
are  multiplied  by  a  factor  of  8.33  the  B  field  is  then  obtained  in  gammas. 
The  electric  field  is  read  directly  in  volts/meter. 

rf'he  incident  slow  wave  disturbance  is  primarily  transverse 
magnetic  at  the  top  of  the  ionosphere  and  becomes  transverse  electric 
in  the  ground.  In  the  region  above  1000  km.  a  standing  wave 
character  is  exhibited  and  is  due  to  reflection  from  file  interface  separating 
the  anisotropic  region  above  1000  km.  from  the  quasi-isotropic  region 
below  1000  km.  Standing  waves  are  also  seen  in  the  region  between 
150  and  1000  km.  and  are  due  to  reflection  within  this  isotropic  region. 

Standing  waves  are  still  discernible  down  through  the  D  region 
and  neutral  atmosphere  but  are  of  considerably  longer  wavelengths  than 
were  seen  in  the  F  region.  The  y  component  of  the  magnetic  field 


1-62 


WAVE 

FIE 


decreases  and  the  z  component  of  the  electric  field  increases  through 
the  D  region  due  to  coupling  caused  by  the  localized  disturbance. 

In  the  neutral  atmosphere  the  z  component  of  the  electric  field  persists 
and  raust  result  in  charge  separation  between  the  lower  ionosphere 
and  ground  level . 

The  x  component  of  the  magnetic  field  grows  gradually  throughout 
the  F  region  and  becomes  nearly  equal  to  the  y  component  in  the  E  region. 
The  x  and  y  raagnetic  field  components  are  again  separated  in  the  D 
region  due  to  coupling.  Inspection  of  the  magnetic  component  change 
with  altitude  reveals  that  whenever  the  x  -  component  is  dominant  on 
the  ground,  the  y  -  component  is  dominant  at  40Cj  km. 

The  electric  field  has  a  similar  behavior  with  the  x  and  y 
components  being  well  separated  in  the  D  and  F  regions  because  of 
coupling.  In  the  E  region  the  x  and  y  electric  field  components  are 
nearly  equal  due  to  the  E  region  anisotropy. 

Within  the  Earth  the  electric  and  magnetic  fields  are  horizontal, 
suggesting  the  properties  of  a  conductor.  The  skin  depth  for  these 
fields  within  the  Earth's  surface  is  less  than  one  km. 

The  wave  fields  for  the  Incident  fast  wave  are  depicted  in  Figure  4.6. 
At  the  top  of  the  ionosphere  the  disturbance  is  neither  transverse  magnetic 
nor  transverse  electric.  However,  the  disturbance  is  transverse  electric 
at  ail  altitudes  below  110  km.  As  was  true  of  the  incident  slow  wave, 
standing  waves  are  discernible  above  250  km.  with  the  magnetic  field 
envelope  reaching  a  maximum  around  the  F  region  ele  tron  density 
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maximum.  The  field  components  become  well  separated  in  the  D  region 
due  to  the  continuation  of  the  TE  mode  in  free  space.  At  ground  level  the 
magnetic  field  is  primarily  in  the  x  direction  and  the  electric  field  .  n 
the  y  direction.  Unlike  the  incident  slow  wave  no  electric  field  7 
component  is  generated  in  the  neutral  atmosphere.  Therefore  the  inciden 
fast  wave  will  not  generate  charge  separation  between  the  ground  and  tlu 
lower  ionosphere.  Also,  the  electric  field  is  attenuated  by  five  orders 
of  magnitude  between  4000  km.  and  the  ground.  Therefore,  the  electric- 
field  of  a  disturbance  generaged  at  4000  km.  may  be  difficult  to  observe 
on  the  ground . 

The  disturbance  within  the  Earth  is  also  transverse  electric. 

4.1S  Wave  Field  Polarization 

The  wave  field  polarizations  for  the  incident  slow  and  fast  waves 
are  given  in  Figures  4.7  and  4.6,  respectively.  The  labeling  of  these 
graphs  is  explained  in  section  3.25. 

T1  e  electric  tiald  of  the  incident  slow  wave  is  left  elliptical!*- 
polarized  in  the  xy  plane  at  4000  km.  with  the  major  axis  in  the  x 
direction.  The  electric  field  of  the  inciden:  fast  wave  is  also  elliptical!* 
polarized  at  4pqq  km.  but  in  the  right  hand  se.  se.  Thus  the  slow 
and  fast  waves  at  4000  km.  may  be  discriminated  by  observation  of  the 
electric  field  polarization.  The  slow  and  fast  modes  at  4000  km.  may 
also  be  differentiated  by  observing  the  magnetic  f*cld  polarization, 
which  is  left  elliptically  polarized  for  the  incident  slow  wave  and  right 
elliptically  polarized  for  the  fast  wave.  However,  or.  the  ground  the 
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Figure  4.7 


QC 


Figure  4 .8 


electric  and  magnetic  field  polarizations  in  the  xy  plane  are 
linear  in  the  same  fashion  for  both  the  incident  slow  and  incident  fast 
waves .  The  two  disturbances  may  be  identified  on  the  ground  only 
by  observation  of  the  z  component  of  the  electric  field,  which  is  zero 
for  the  incident  fast  wave  and  non-zero  for  the  incident  slow  wave. 

The  electric  and  magnetic  fields  of  both  the  incident  slow  and 
fast  waves  undergo  a  polarization  reversal  in  the  xy  plane  at 
approximately  1000  km.  This  polarization  reversal  is  due  to  the  sudden 
change  of  the  medium  from  anisotropic  above  1000  km.  to  quasi-isotropic 
below  1000  km.  Also  the  magnetic  field  of  the  incident  slow  wave  changes 
back  to  left  elliptically  polarized  in  the  xy  plane  at.  90  km. ,  which 
corresponds  to  the  point  at  which  the  x  component  of  the  magnetic  field 
becomes  larger  than  the  y  component. 

The  incident  fast  wave  electric  field  is  left  circularly  polarized 
between  1000  km.  and  100  km.  and  then  becomes  linearly  polarized  in 
the  y  direction.  The  incident  slow  wave  electric  field  is  right 
elliptically  polarized  between  1000  km.  and  ground  surface,  at 
which  it  becomes  linearly  polarized  in  the  y  direction  in  the  xy  plane. 

Experimental  observations  indicate  that  at  high  latitudes  the  Pc  1 
signals  received  on  the  ground  during  the  day  are  primarily  left 
elliptically  polarized  [Campbell,  1965,  1967,  1968]  but  that  frequently 
the  signals  are  right,  linearly,  or  irregularly  polarized  [Heacock,  1970], 
Some  evidence  [Heacock,  1970]  suggests  that  Pc  1  events  are  a  super¬ 
position  of  two  or  more  events  and  Pope  [1964]  has  attempted  to  explain 
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the  polarization  behavior  of  Pc  1  events  in  terms  of  a  superposition 
of  successive  rising  tone  signals,  the  resultant  polarization  sense  of 
which  is  determined  by  the  relative  magnitudes  of  the  signals  and  the 
relative  inclination  of  their  polarization  axes.  More  recently  Campbell 
[1967]  has  concluded  from  observational  evidence  of  Pc  1  attenuation 
characteristics  and  excellent  geomagnetic  conjugacy  that  the  signals 
received  at  ground  level  are  L  waves  having  R  wave  characteristic 
absorption. 

Experimental  workers  (e.g.,  Heacock  [1970])  indicate  a 
complexity  in  the  polarization  of  Pc  1  signals  observ  ed  at  auroral  zone 
stations.  The  results  of  the  present  work  indicate  that  the  first 
independent  wave,  having  left  handed  elliptical  polarization  at  the  top 
of  the  ionosphere,  shows  less  attenuation  in  the  course  of  ionospheric 
transmission  than  the  second  independent  wave,  the  polarization  of 
which  is  right  handed  and  almost  circular  at  the  top  oi  .he  ionosphere. 
This  result  agrees  with  Campbell's  implication,  which  was  in 
contradiction  with  the  expectation  from  the  case  of  an  infinite  plane 
wave  in  a  uniform  plasma.  Although  this  agreement  comes  essentially 
from  the  consequence  of  magneto-ionic  coupling  the  actual  process  is 
not  so  simple. 

In  the  present  model  of  the  ionosphere,  together  with  the 
localization  of  the  disturbance  in  the  x-direction,  mode  interchange 
happens  at  slightly  above  1000  km.  and  at  150  km.  Between  these 
levels,  the  medium  is  essentially  isotropic,  so  that  there  is  no 
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appreciable  difference  in  the  attenuation  of  two  characteristic  waves. 


The  differentiation  of  attenuation  occurs  only  in  the  E  and  D  regions . 

As  seen  from  Figure  4.4  the  incident  "L"  wave  consists  mainly  of  the 
slow  mode,  whereas  the  incident  "R"  wave  is  comprised  of  both  the 
fast  and  slow  modes .  By  virtue  of  the  mode  coupling  around  the  upper 
mode-interchange  level,  the  incident  "L"  wave  becomes  a  composite 
wave  of  slow  and  fast  modes  with  the  same  magnitude  in  the  upper 
F  region,  while  the  incident  "R"  wave  becomes  a  rather  pure  slow  mode 
in  that  region.  Then,  a  mode  interchange  again  occurs  at  150  km. 

Below  this  level,  the  differentiation  in  the  attenuation  of  different 
modes  becomes  obvious.  As  a  result,  the  incident  "R"  wave  is  more 
attenuated,  compared  to  the  incident  "L"  wave. 

In  contrast  with  Heacock's  observational  findings,  which  show  a 
change  of  the  rotation  sense  of  wave  polarization  from  one  location  to 
another  and  from  one  occasion  to  another,  the  present  analysis  gives 
only  linear  polarization  in  the  x-direction  on  the  ground.  The  same 
result  was  obtained  at  different  frequencies  in  the  ULF  range.  This  is 
a  consequence  of  the  assumption  that  the  disturbance  is  localized  only 
in  the  x-direction.  In  other  words  the  one  dimensional  localization  is 
equivalent  to  the  specification  of  the  dominant  TE  mode  in  the  free 
space  between  the  ionosphere  and  the  Earth's  surface.  From  this 

analysis  the  conjecture  may  be  made  that  if  the  disturbance  function 

-Xx  ifcjmu 

is  of  the  form  C  •  £  a  *  the  resultant  polarization  can  be  elintical 
and  the  sense  of  rotation  is  determined  by  the  sign  of  the  ratio  Va  • 
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Also,  there  is  observational  evidence  showing  the  structural 
ionization  across  the  auroral  zone,  in  particular,  during  auroral  events, 
so  that  Pc  1  wave  disturbances  may  be  confined  in  the  vicinity  of  the 
structural  ionization,  such  as  the  auroral  column  and/or  the  auroral 
arc.  On  some  occasions,  Pc  1  waves  can  be  trapped  between  two 
sheets  of  the  auroral  arcs.  Therefore,  in  order  to  clarify  the  spatial 
change  of  the  polarization  rotation  sense,  the  wave  guide  mode 
propagating  along  the  structural  ionization  observed  in  the  auroral  zone 
must  be  examined  by  taking  into  account  a  certain  bi-directional 
disturbance  function.  Thus,  as  far  as  auroral  zone  phenomena  are 
concerned,  the  assumption  of  an  infinite  plane  wave  in  a  uniform  plasma 
is  not  sufficient  to  account  for  various  complexities  in  the  nature  of 
the  phenomena. 

Also,  the  present  results  give  larger  attenuation  for  both  the 
incident  "L"  and  "R"  waves  than  given  by  previous  investigations. 

In  addition,  the  present  work  shows  a  greater  difference  between  the 
incident  "L"  and  "R"  wave  attenuations  than  predicted  by  previous 
investigators.  These  differences  are  attributed  to  the  inclusion  of  the 
ion  resonant  charge  exchange  processes. 

Theoretical  studies  [Altman  and  Fijalkow,  1963,  1969;  Altman  and 
Cory,  196^]  have  suggested  that  intermode  coupling  channels  energy 
from  L  to  R  modes  in  the  upper  E  region.  These  studies  have  concluded 
that  the  coupling  increases  as  the  angle  between  the  wave  normal  and  the 
magnetic  field  increases.  Heacock  [1970]  reconciles  his  dual-site 
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(College  and  Kotzebue,  Alaska)  polarization  differences  by  invoking 
Altman  and  Fijalkow's  L-to-R  coupling,  assuming  that  Colleye  and 
Kotzebue  are  both  at  the  feet  of  the  disturbance  source  lines,  and 
assuming  that  flaring  out  of  the  L-wave  ray  bundle  occurs  as  the 
waves  propagate  through  the  F  region.  Large  amplitude  L  sense 
should  then  be  observed  near  the  center  of  the  disturbed  field  line 
region  and  the  polarization  should  change  gradually  to  R  sense  with 
decreasing  amplitude  as  the  distance  increases  radially  from  the 
center  of  the  disturbed  region. 

4.2  Localized  Disturbance  at  1/3  c/s. 

4.21  Booker  Quartic  Roots 

Table  4.2  summarizes  the  properties  of  the  slow  and  fast 
wave  disturbances  at  1/3  c/s  in  the  various  regions  of  the  Earth's 
environment.  The  Booker  quartic  roots  are  presented  in  Figure  4.9. 
The  labeling  of  these  cuivos  is  the  same  as  that  used  in  section 
3.21.  The  actual  values  of  the  indices  of  refraction  for  propagation 
in  the  ^-direction  are  obtained  by  multiplying  the  values  obtained 
from  the  curves  by  a  factor  of  22.5. 

Above  600  km .  the  region  is  anisotropic  and  the  slow  wave  is 
propagating  and  the  fast  wave  is  evanescent.  From  magneto-ionic 
theory  one  would  expect  that  both  the  slow  and  fast  modes  would  be 
propagating  through  this  region  because  the  wave  frequency  is  less 
than  the  ion  cyclotron  gyrofrequency .  However,  the  dispersion 
relation  given  in  section  2.2  shows  that  above  600  km.  the  Booker 
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TABLE  4.2  Medium  characteristics  at  1/3  c/s  for  a  localized 
disturbance  incident  at  an  oblique  angle  from  the  vertical  at  5000  Km. 
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quartlc  roots  for  the  fast  and  slow  waves  are  given  approximately  by 


respectively,  where X is  the  Fourier  number. 


S"  = 


,  and  r  = 


XIV  ftt 

~rr 


At  500  km.  these  waves  are  given  by 

Vt^l? 


Thus  the  imaginary  part  of  the  Pedersen  conductivity  increases 

P2 

sufficiently  to  overcome  the  localization  factor**-  and  enables  the 
fast  wave  to  become  propagating .  This  effect  holds  true  for  the 
region  200-600  km.,  which  is  quasi-isotropic  and  both  modes  are 
propagating.  The  real  parts  of  the  quartic  roots  become  nearly 
equal  to  the  imaginary  parts  and  this  region  therefore  is  similar  to 
a  metal .  Later  in  the  analysis  this  conducting  region  will  be  seen 
to  serve  as  a  waveguide  for  rhe  fast  mode. 

Between  about  200  km.  and  70  km.  the  slow  mode  is  propagating 
and  the  fast  mode  is  evanescent.  Therefore  the  propagating  region 
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for  the  fast  mode  is  sandwiched  between  two  strongly  evanescent 
regions  and  energy  trapping  of  the  fast  wave  occurs  between  approxi¬ 
mately  200  and  600  km . 

Charge  exchange  has  basically  the  same  effect  on  the  quartic 
roots  at  0.333  c/s  as  at  3.0  c/s.  However,  at  the  lower  frequency 
the  coupling  coefficients  are  becoming  comparable  in  magnitude 
to  the  quartic  roots  and  may  somewhat  mask  the  effects  of  charge 
exchange.  Nevertheless  some  absorption  still  takes  place  between 
about  160  km.  and  220  km.  where  the  quartic  roots  are  still  larger 
than  the  coupling  coefficients. 

Both  the  slow  and  fast  modes  are  evanescent  in  the  quasi¬ 
isotropic  portion  of  the  D  region  below  70  km.  The  neutral  atmosphere 
and  the  Earth's  interior  are  Isotropic.  The  quartic  roots  for  the  neutral 
atmosphere  are  once  again  pure  imaginary  due  to  the  localized  nature  of 
the  disturbance. 

4.22  Magneto-ionic  Coupling  Coefficients 

Figure  4. 10  shows  the  altitude  dependence  of  the  magneto- ionic 
coupling  coefficients  for  the  localized  disturbance  at  1/3  c/s.  These 
curves  have  the  same  labeling  as  those  analyzed  in  section  3.22. 
Reflection  type  coupling  is  most  important  down  to  300  km.  Coefficients 
for  the  reflection  of  the  slow  wave  are  greater  than  those  for  the 
reflection  of  the  fast  wave  above  about  700  km .  The  reflection 
coefficients  for  the  fast  mode  increase  rapidly  l>etween  1000  km.  and 
600  km.  due  to  the  rapid  change  in  the  medusr  character 


1-77 


MAGNETO- IONIC  COUPLING 


from  highly  anisotropic  to  isotropic  in  the  region  between  200 
and  600  km.  The  coupling  coefficients  for  reflection  of  the  slow 
wave  are  the  same  order  of  magnitude  as  the  slow  wave  quartic 
roots  throughout  the  F  region.  The  coupling  coefficients  for 
reflection  of  the  fast  wave  are  smaller  than  the  fast  wave  quartic  roots 
above  700  km.  but  become  larger  than  the  quartic  roots  in  the  region 
200-600  km.  Thus  energy  from  the  fast  wave  is  trapped  in  this 
conducting  region,  which  behaves  as  a  waveguide  for  the  fast  wave. 
Mode  coupling  continuously  increases  throughout  the  F  region  and 
reaches  a  maximum  around  the  F  region  maximum.  Mode  reflection  type 
coupling  also  appears  in  the  F  region  but  is  of  less  importance  than  pure 
mode  coupling.  Once  again  in  the  E  region  reflection  type  coupling 
is  important.  In  fact,  for  both  modes  the  reflection  coupling  coefficients 
are  larger  than  the  quartic  roots.  There  is  little  mode  coupling  in  the 
E  region  as  evidenced  by  the  good  separation  of  the  quartic  roots  for  the 
siow  and  fast  modes.  Mode  coupling  and  mode  reflection  type 
coupling  both  become  much  larger  than  the  quartic  roots  in  the  lower 
D  region.  Thus  mode  coupling  is  important  in  the  D  and  F  regions  and 
unimportant  in  the  E  region.  This  result  seems  to  be  in  contradiction 
to  the  suggestion  made  by  Heacock  [1970] . 

Reflection  type  coupling  occurs  with...  the  surface  layer  of  the 
Earth  and  is  a  maximum  at  a  depth  of  about  2.5  km.  This  reflection 
region  corresponds  to  the  rapid  change  in  the  earth  conductivity.  In 
this  region  the  coupling  coefficients  are  smaller  than  the  quartic  roots 
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but  are  of  comparable  magnitude  at  a  depth  of  2.5  km. 

4.23  Characteristic  Wave  Mixing  Ratios 
The  mixing  ratios  for  the  incident  slow  and  fast  waves  at 
1/3  c/s  are  given  in  Figure  4.11.  The  labeling  of  these  curves  is  the 
same  as  that  used  in  section  3.23. 

For  the  incident  slow  wave  the  downgoing  and  upgoing  slow  waves 
are  of  similar  magnitude  throughout  the  F  region,  indicating  that 
reflection  is  continuously  taking  place.  The  downgoing  and  upgoing 
fast  waves  are  gradually  generated  by  mode  coupling .  In  the  D  and 
E  regions  the  upgoing  and  downgoing  waves  are  of  similar  magnitude 
which  indicates  that  standing  waves  are  formed  and  that  the  D  and  E 
regions  behave  as  a  waveguide  for  the  incident  slow  wave.  In  the 
neutral  atmosphere  the  transverse  electric  mode  is  dominant  over 
the  transverse  magnetic  mode.  Also  almost  complete  reflection  takes 
place  at  ground  surface.  Thus  the  Earth  is  a  good  conductor  even 
at  such  a  low  frequency.  Inside  the  Earth  the  downgoing  transverse 
electric  mode  is  dominant  but  the  upgoing  transverse  electric  mode 
is  also  present  due  to  reflection. 

The  mixing  ratios  for  the  fast  mode  decrease  by  three  orders  of 
magnitude  between  1000  km.  and  1000  km.  and  by  four  orders  of 
magnitude  oetween  4000  km.  and  ground  level.  Therefore  even  if  a 
disturoance  arrives  as  a  fast  wave  at  4000  km.  it  might  not  be  seen  on 
the  ground.  If  the  fast  mode  is  generated  by  a  disturbance  below 
1000  km.  the  wave  is  more  likely  to  be  observed  on  the  ground.  Below 
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600  km.  the  downgoing  and  upgoiny  waves  become  similar  in  magnitude 
and  show  that  the  region  between  about  250  km.  and  600  km.  serves  as 
a  waveguide  for  the  fast  wave.  In  the  neutral  atmosphere  the  transverse 
electric  mode  is  the  only  mode  present  and  almost  complete  reflection 
of  this  mode  takes  place  on  the  ground.  The  downgoing  transverse 
electric  mode  is  dominant  in  the  Earth  but  the  upgoing  transverse 
electric  mode  is  also  present  to  a  small  degree  due  to  reflection. 

4.24  Wave  Fields 

The  wave  fields  of  the  incident  slow  mode  are  shown  in 
Figure  4.12  .  The  magnitudes  of  the  electric  and  magnetic  fields  are 
all  relative  to  the  values  at  4000  km.  If  the  magnetic  field  values  are 
multiplied  by  a  factor  of  75  the  B  field  is  then  obtained  in  gammas.  The 
electric  field  is  read  directly  in  volts/meter. 

The  incident  slow  mode  disturbance  is  transverse  magnetic  at 
the  top  of  the  ionosphere  and  becomes  transverse  electric  within  the 
ground.  The  y  component  of  the  magnetic  field  increases  through  the 
F  region  and  reaches  a  maximum  around  300  km.  Magnetic  energy  is  being 
stored  in  the  region  200-600  km.  The  y  component  of  the  magnetic 

4 

field  then  decreases  in  may*.  j  through  the  D  and  E  regions .  The 
x  and  z  components  of  the  magnetic  field  increase  through  the  F  region 
due  to  mode  coupling  and  exhibit  somewhat  of  a  standing  wave  character 
below  about  300  km .  The  x  and  z  components  of  the  magnetic  field 
become  dominant  over  the  y  component  below  100  km.  The  x  component 
of  the  magnetic  field  is  seen  in  Figure  4.12  to  be  dominant  on  the 


1-82 


-S  -J  -l  0  1  -1  -t  -J  -«  -J 

H  FIU.0  CO«W»tJ  E  MELO  CQPrONfNIS 

8««.  3  »oo( 


Figure  4.12 


ground  whereas  the  y  component  is  dominant  at  4000  km. 

The  electric  field  is  in  the  x  direction  at  the  top  of  the  ionosphere 
but  in  the  z  direction  at  ground  level.  The  y  component,  grows  in 
magnitude  through  the  F  region  and  becomes  equal  to  the  x  component 
around  150  km.  Below  this  level  the  x  and  y  components  exhibit  a 
standing  wave  formation.  Below  about  90  km.  the  x  and  y  components 
become  somewhat  separated  in  the  D  region.  The  z  component 
also  increases  in  magnitude  through  this  low  conducting  region  and 
becomes  dominant  in  the  neutral  atmosphere.  Therefore  charge  separation 
must  occur  between  the  Earth’s  surface  and  the  lower  ionosphere.  Both 
the  electric  and  magnetic  fields  decrease  logarithmically  within  the 
Earth.  Even  at  this  low  frequency  the  Earth  is  a  good  conductor. 

The  wave  fields  for  the  incident  fast  wave  are  shown  in  Figure 
4.13.  Both  the  electric  and  magnetic  fields  of  this  disturbance  are 
attenuated  by  three  orders  of  magnitude  between  4000  and  1000  km. 

The  magnetic  field  is  attenuated  by  four  orders  of  magnitude  between 
4000  km.  and  ground  level  and  the  electric  field  is  attenuated  by  seven 
orders  of  magnitude  in  the  same  altitude  range .  Therefore  an  incident 
fast  wave  arriving  at  4000  km.  may  be  difficult  to  observe  on  the 
ground.  A  fast  mode  disturbance  generated  below  1000  km.  would  be 
more  easily  observed.  The  incident  fast  wave  is  transverse  electric 
throughout  the  ionosphere,  neutral  atmosphere  and  the  Earth's  interior, 
with  the  exception  of  the  region  between  200  and  600  km.  altitude 
where  strong  coupling  causes  the  x  and  y  components  of  the  electric 
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Figure  4.13 


field  to  be  nearly  equal.  A  z  component  of  the  electric  field  is 
generated  in  the  lower  D  region  but  unlike  the  incident  slow  mode 
the  z  component  of  the  electric  field  does  not  appear  in  the  neutral 
atmosphere  or  on  the  ground. 

4.25  Wave  Field  Polarization 

The  wave  field  polarizations  for  the  incident  slow  and  fast 
waves  are  given  in  Figures  4.14  and  4.15,  respectively.  These 
figures  are  interpreted  in  the  same  manner  outlined  in  section  3.25. 

The  polarization  characteristics  of  the  resultant  waves  at  this  frequency 
(0.333  c/s),  which  is  near  the  lower  limit  of  the  Pc  1  frequency  range, 
are  also  in  agreement  with  Campbell's  [1967]  suggestion  that  the  signal 
observed  on  the  ground  represents  a  disturbance  which  has  propagated 
through  the  ionosphere  with  R  mode  attenuation.  The  incident  slow 
mode  magnetic  field  is  linearly  polarized  down  to  250  km.  where  it 
then  becomes  right  ellipticaliy  polarized.  At  110  km.  the  magnetic 
field  changes  to  left  ellipticaliy  polarized  and  then  to  linearly  polarized 
in  the  D  region.  The  electric  wave  field  of  the  incident  slow  mode  it 
linearly  polarized  down  to  250  km.  and  then  becomes  right  ellipticaliy 
polarized,  remaining  as  such  down  to  ground  level  where  the  polarization 
is  linear.  Thus  the  incident  slow  mode  basically  has  R  mode  charac¬ 
teristics  and  is  less  attenuated  than  the  incident  fast  wave. 

The  incident  fast  wave  magnetic  field  is  linearly  polarized  at 
4000  km. ,  changes  to  left  ellipticaliy  polarized  at  1000  km.,  and 
remains  left  ellipticaliy  polarized  down  to  the  D region  where  it  becomes 
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Figure  4.14 


Figure  4.15 


linearly  polarized  again.  The  wave  electric  field  has  the  same 
polarization  characteristics  as  the  magnetic  field  in  the  xy  plane 
except  that  the  major  axes  of  the  polarization  curves  are  approximately 
at  90°  to  each  other.  Therefore  the  incident  fast  mode  is  found  to  be 
in  the  L  mode  throughout  most  of  the  ionosphere  and  is  more  attenuated 
than  the  incident  slow  wave . 

If  the  bi-directional  disturbance  function  is  assumed  to  be  of 
the  form  £  t  an  elliptically  polarized  wave  may  then 

appear  at  ground  level . 

4.3  Localized  Disturbance  at  1/300  c/s 

4.31  Booker  Quartic  Roots 

Table  4.3  summarizes  the  properties  of  the  slow  and  fast  wave 
localized  disturbances  at  1/300  c/s  in  the  various  regions  of  the  Earth’s 
environment.  Figure  4.16  shows  the  altitude  dependence  of  the  Booker 
quartic  roots .  The  indices  of  refraction  for  propagation  in  the  z-directlon 
are  obtained  by  multiplying  the  values  obtained  from  the  curves  in 

3 

Figure  4.16  by  a  factor  of  2.25  x  10  . 

The  quartic  roots  reveal  that  the  ionosphere  Is  anisotropic 
throughout.  The  slow  wave  is  propagating  down  to  about  80  km.  where 
it  becomes  evanescent.  The  fast  wave  is  highly  evanescent  throughout 
the  entire  Ionosphere.  This  is  similar  to  the  fast  wave  disturbance  at 
1/3  c/s,  which  was  evanescent  above  600  km.  The  same  reason  given 
to  exDlrin  the  evanescence  of  the  fast  wave  in  the  Pc  1  frequency  range 
may  also  be  given  to  explain  the  evanescence  of  the  fast  wave  in  the 
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TABLE  4.3  Medium  characteristics  at  1/300  c/s  for  a  localized 
disturbance  incident  at  an  oblique  angle  from  the  vertical  at  S000  Km. 


Pc  5  frequency  range.  The  dispersion  relation  given  in  section  2 .2 


predicts  that  above  2C0  km.  the  3ooker  quartic  roots  for  the  fast 
and  slow  waves  are  given  approximately  by 

%  -  ‘I 
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Thus  the  disturbance  factor  jfc  is  sufficiently  large  to  cause  the  fast 
wave  to  become  evanescent. 

At  approximately  200  km.  the  real  and  imaginary  parts  of  the 
quartic  roots  for  the  slow  wave  become  nearly  equal .  Also  the  real 
part  of  the  quartic  root  of  the  fast  wave  increases  rapidly,  reaching 
a  maximum  at  about  150  km.  but  remaining  a  factor  of  two  less  than 
the  imaginary  part.  Consequently  mode  coupling  becomes  large  in 
this  region. 

From  the  dispersion  relation  it  may  be  shown  that  at  200  km. 
the  fast  and  slow  waves  are  approximated  by 

The  rapid  increase  of  the  real  part  of  the  Pedersen  conductivity  with 
decreasing  altitude  causes  the  real  part  of  the  fast  wave  quartic  root 

to  approach  the  imaginary  part,  but  the  increase  in  ^  is  not 

IX 
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sufficient  to  overcome  the  disturbance  factor^  and  to  allow  the  fast 
wave  to  become  propagating. 

Between  150  km.  and  600  km.  is  enhanced  by  as  much 

.*  T 

as  an  order  of  magnitude  when  the  effect  of  charge  exchange  of 
G+with  Ois  included.  Therefore  the  effect  of  charge  exchange  is  co 
cause  the  slow  wave  to  become  more  propagating  and  the  fast  wave  to 
approach  more  rapidly  toward  propagating  conditions . 

The  real  and  imaginary  parts  of  the  quartic  roots  for  the  slow 
wave  again  become  unequal  in  the  E  region  below  150  km.  so  that  the  slow 
wave  Is  highly  propagating  in  this  region. 

Once  again  the  quartic  roots  are  pure  imaginary  ir.  the  neutral 
atmosphere  due  to  the  presence  of  localization  in  the  horizontal  direction 
of  the  disturbance.  Also,  the  interior  of  the  Earth  is  like  a  metal. 

4.32  Magneto-ionic  Coupling  Coefficients 

Figure  4.17  shows  the  altitude  dependence  of  the  magneto-ionic 
coupling  coefficients  for  the  localized  disturbance  at  1/300  c/s.  These 
curves  have  the  same  labeling  as  those  analyzed  in  section  3.22. 

Reflection  type  coupling  for  both  the  slow  and  fast  waves  is  seen 
to  be  the  dominant  type  of  coupling  above  the  F  region.  In  fact 
the  coupling  coefficients  are  everywhere  larger  than  Re  [qs] , 
indicating  that  the  coupling  coefficients  are  more  important  than  the 
quartic  roots  in  determining  the  final  wave  characteristics .  In  the  region 
150-250  km.  mode  coupling  and  mode  reflection  type  coupling  become 
important  due  to  the  fact  that  the  medium  is  becoming  isotropic. 
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As  will  be  shown  later,  across  this  mode  coupling  region  the  incident 
slow  wave  changes  from  quasi-TM  at  the  higher  altitudes  to  quasi-TE 
at  the  lower  altitudes. 

In  the  E  region  reflection  type  coupling  once  again  dominates  for 
the  slow  wave.  Mode  coupling  becomes  very  large  in  the  lower  D 
region,  which  was  also  true  at  higher  frequencies  and  was  shown  to  be 
caused  by  the  localization  of  the  disturbance.  For  the  fast  wave  mode 
coupling  is  most  important  throughout  the  lower  isonosphere. 

Reflection  type  coupling  occurs  within  the  surface  layer  of  the 
Earth  between  1.2  km.  and  3.5  km.  depth  and  is  a  maximum  at  a  depth 
of  2.5  km.  At  all  depths  within  this  reflection  layer  the  reflection 
type  coupling  coefficients  are  larger  than  the  quartic  roots .  This 
reflection  region  is  caused  by  the  rapid  change  in  the  earth  conductivity. 
The  mixing  ratios  reveal  that  a  reflected  TE  wave  is  generated  by  the 
reflection  coupling . 

4.33  Characteristic  Wave  Mixing  Ratios 

The  characteristic  wave  mixing  ratios  for  the  incident  slow  and 
fast  waves  arc  given  in  Figure  4.18.  These  curves  have  the  same 
labeling  as  those  analyzed  in  section  3.23. 

For  the  incident  slow  wave  the  downgoing  and  upgoing  slow  waves 
are  approximately  equal  above  the  F  region,  indicating  that  reflection  is 
continuously  taking  place.  The  downgoing  and  upgoing  fast  waves  are 
generated  around  the  F  region  max  as  a  result  of  mode  and  mode  reflection 
type  coupling.  After  passing  through  the  strong  mode  coupling  region, 
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Figure  4.18 


all  the  characteristic  waves  become  the  same  order  of  magnitude .  Below 
150  km.  in  the  D  and  E  regions  mixing  of  the  characteristic  waves  is 
quite  good  due  to  mode  coupling  and  mode  reflection  type  coupling. 

In  the  neutral  atmosphere  the  waves  are  mostly  transverse  magnetic 
but  with  a  small  amount  of  transverse  electric  waves  present.  The 
surface  of  the  Earth  acts  as  a  good  conductor  even  at  Pc  5  frequencies , 
although  the  fields  are  able  to  penetrate  more  deeply  than  at  the  higher 
frequencies.  Within  the  Earth  the  resultant  wave  is  primarily  downgoing 
transverse  electric . 

The  incident  fast  wave  consists  primarily  of  the  downgoing  fast 
characteristic  wave  at  the  top  of  th  ionosphere.  This  character 
predominates  down  to  below  the  F  region  maximum,  because  of  the 
constancy  of  Im  [q^]  throughout  the  higher  region.  Around  the  F  region 
maximum  the  slow  characteristic  waves  are  generated  by  the  mode  and 
mode  reflection  type  coupling.  At  higher  altitudes  a  similar  behavior 
is  seen  in  the  variation  of  the  slow  wave  components  contained  in  the 
initial  slow  and  fast  modes.  Therefore,  the  slow  wave  components  in 
the  incident  fast  wave  can  be  thought,  of  as  being  included  at  incidence. 
The  upgoing  fast  wave  originates  around  the  F  region  maximum  where 
reflection  type  coupling  is  strong.  The  incident  fast  wave  is  highly 
attenuated  and  has  decreased  in  magnitude  by  three  orders  of  magnitude 
at  1000  km. 

Although  the  mixing  ratios  in  the  neutial  atmosphere  show  a 
predominance  of  the  TM  mode,  the  reflection  of  the  TM  mode  at  ground 
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surface  is  almost  perfect  in  the  sense  that  the  downgoing  and  upgoing 
components  have  the  same  magnitude  and  are  180°  out  of  phase.  As  a 
result,  the  TM  upgoing  and  downgoing  waves  cancel  each  other.  On 
the  other  hand,  while  the  TE  mode  mixing  ratios  are  much  smaller  than 
those  for  the  TM  mode,  the  reflection  of  the  TE  mode  at  ground  surface 
is  not  perfect;  the  incoming  and  reflected  waves  are  almost  in  phase 
and  additive.  Consequently,  the  TE  upgoing  and  downgoing  waves 
interfere  constructively  in  free  space,  instead  of  destructively,  as  is 
the  case  of  the  TM  mode.  Therefore,  the  TE  -  like  wave  field  remains 
in  the  neutral  atmosphere,  as  will  be  seen  later. 

Within  the  ground  the  resultant  wave  is  transverse  electric  with  some 
reflection  taking  place.  This  mode  also  penetrates  to  greater  depths  at 
this  frequency  than  at  higher  frequencies. 

4.34  Wave  fields 

The  wave  fields  of  the  incident  fast  wave  are  shown  in  Figure  4.19. 
The  magnitudes  of  the  electric  and  magnetic  fields  are  all  relative  to 
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the  values  at  4000  km.  If  the  magnetic  iield  values  are  multiplied  by  a 
factor  of  7500  the  B  field  is  then  obtained  in  gammas.  The  electric 
field  is  read  directly  in  volts/meter. 

The  incident  fast  wave  is  essentially  transverse  electric  throughout 
the  entire  region  of  consideration.  No  obvious  standing  waves  are 
visible  in  the  wave  field  envelopes  as  there  were  at  the  higher  frequencies 
considered  in  this  investigation.  The  wave  fields  are  already  attenuated 
by  three  orders  of  magnitude  at  1000  km.  At  ground  level  the  magnetic 
and  electric  fields  are  attenuated  by  four  and  five  orders  of  magnitude, 
respectively,  compared  to  the  values  at  4000  km.  Consequently  it  seems 
likely  that  the  incident  fast  wave  may  be  difficult  to  observe  at  ground 
level.  A  small  amount  of  electric  field  z-component  is  generated  in  the 
lower  D  region  due  to  the  localized  disturbance  effect,  but  no  rlectrc 
field  z-component  is  generated  in  the  neutral  atmosphere.  Therefore  the 
incident  fast  wave  will  not  produce  charge  separation  between  the  lower 
ionosphere  and  ground  level. 

The  presence  of  the  z  component  of  the  magnetic  field  throughout 
the  entire  region  of  consideration  is  also  a  consequence  of  the  localization 
of  the  disturbance.  In  other  words  this  mode  is  mostly  of  Tf,‘  character  and 
the  magnetic  z-component  persists  everywhere.  For  the  plane  wave  case  a 
z-component  of  the  magnetic  field  does  not  appear. 

1.;  the  ground  the  fields  are  attenuated  logarithmically  but  penetrate 
to  greater  depths  than  at  higher  frequencies . 

The  wave  fields  of  the  incident  slow  wave  are  depicted  in  Figure 
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4  .20.  This  disturbance  is  transverse  magnetic  above  the  F  region 
maximum  and  is  transformer  into  transverse  electric  below  the  ionosphere. 
The  transition  takes  place  between  the  F  region  maximum  and  100  km. 
Within  the  transition  layer  the  magnetic  field  character  changes  from 
TM  to  TE  due  to  the  presence  of  mode  coupling.  However,  the  character 
of  the  electric  field  does  not  change  across  the  transition  layer.  The 
electric  field  mode  change  occurs  abruptly  at  the  lower  boundary  of 
the  ionosphere. 

Within  the  transition  layer  between  200  km.  and  100  km.  a 
horizontal  electric  current  flows  rather  strongly  and  gives  rise  to  a 
steep  change  in  the  x  and  y  magnetic  components.  This  electric  current 
iayer  is  necessary  for  the  transition  from  TM-like  at  higher  altitudes  to 
TE-like  at  lower  altitudes.  Also  it  should  be  stressed  here  that  whenever 
the  disturbance  is  specified  by  a  large  localizing  factor  in  the  x-direction, 
the  wave  field  ought  to  behave  as  a  TE  mode  in  the  neutral  atmosphere, 
which  is  bounded  by  the  conducting  earth  and  the  ionosphere.  Moreover, 
it  is  interesting  that  the  electric  current  layer  coincides  with  the 
maximum  conduction  region  of  the  ionosphere  (See  Part  11). 

In  the  lower  D  region,  where  the  conductivity  is  rather  low,  a 
^-component  of  the  electric  field  appears  as  a  result  of  the  localized 
nature  of  the  disturbance.  Th»  electric  field  z-componcnt  persist.- 
throughout  the  neutral  atmosphere  and  must  cause  charge  separation  to 
occur  between  the  lower  ionosphere  and  ground  surface.  This  effect 
is  not  observed  for  the  incident  fast  wave  or  the  incident  plane  wave. 
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Figure  4.20 


The  incident  slow  wave  is  not  severely  attenuated  between  4000  km. 


and  ground  level  whereas  the  incident  fast  wave  is  severly 
attenuated.  In  fact  ..he  value  of  the  electric  field  z-component  of  the 
incident  slow  wave  at  ground  surface  is  enhanced  by  a  factor  of 
four  over  the  x-component  of  the  electric  field  at  4000  km. 

4.35  Wave  Field  Polarization 

In  Figure  4.21  the  incident  fast  wave  is  seen  to  have  the 
magnetic  and  electric  fields  linearly  polarized  in  the  horizontal 
plane  in  the  x  ~nd  y  directions,  respectively,  throughout  the 
ionosphere  and  at  ground  level.  However,  for  the  magnetic  field 
th  xz  plane  is  i  ost  important  because  the  x  and  z  magnetic  field 
components  are  dominant  thraugbout  the  entire  region  of  consideration. 
In  the  xz  plane  the  incident  lust  wave  is  right  circularly  polarized 
at  4000  km.  and  right  elliptically  polarized  at  ground  level.  The 
intervening  ionosphere  serves  to  modify  the  wave  from  circular 
to  elliptical  polarization. 

The  incident  fast  wave  is  seen  to  be  compressive  in  the 
presence  of  the  £  X  S  drift  at  high  altitudes.  As  was 
seen  in  the  wave  fields  in  Figure  4.19,  this  localized  compressive 
mode  is  highly  attenuated  in  the  form  of  an  evanescent  wave,  even 
at  high  altitudes.  Therefore,  the  conclusion  is  drawn  that  a 
localized  disturbance  of  compressive  nature  cannot  be  transmitted 
along  a  magnetic  field  line  down  to  ground  surface  as  long  as  the 
cold  plasma  approximation  is  made. 


1-103 


At  4000  km.  the  wave  normal  direction  for  the  incident  slow  wave 
is  almost  perpendicular  to  the  magnetic  field  line  because  of  the  smallness 
of  Re  [qg] .  At  the  higher  altitudes  this  wave  is  torsional,  i.  e.,  an 
Alfven  wave,  and  is  almost  linearly  polarized  (Figure  4.22).  This 
torsional  localized  disturbance  is  quite  well  transmitted  down  to 
ground  level.  However,  as  already  seen  in  the  wave  fields  (Figure  4.20), 
a  mode  change  from  torsional  at  higher  altitudes  to  rather  compressional 
at  lower  altituces  takes  place  across  the  good  conducting  region  centered 
on  150  km.  Therefore,  the  ionospheric  conducting  region  acts  as  a 
surface  current  layer  separating  the  upper  ionosphere  of  hydromagnetic 
character  from  the  lower  poorly  conducting  ionsophere.  Because  of  the 
strong  localization  in  the  horizontal  direction  and  of  the  difference  in 
the  reflection  at  the  ground  surface  oetween  the  TE  and  TM  modes,  the 
wave  in  free  space  is  forced  to  be  of  TE  character.  The  mode  conversion 
which  takes  place  between  higher  and  lower  regions  gives  rise  to  the 
surface  current  layer  in  the  good  conducting  region  of  the  ionosphere. 

The  incident  fast  wave  maintains  the  TE  character  throughout  the 
entire  region  under  consideration,  so  that  any  obvious  surface  current 
layer  does  not  appear. 

Thus,  from  the  present  analysis,  it  may  be  concluded  that  a 
torsional  disturbance,  even  if  localized,  can  penetrate  from  the 
magnetosphere  down  to  the  ground  surface.  However,  a  localized 
compressional  disturbance  can  not  be  observed  of  the  ground  because 
of  its  evanescent  character  in  the  high  i  icre  and  in  the  magnetosphere. 
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Figure  4.22 


Localized  torsional  waves  resulting  from  the  Kelvin- Helmholtz 
instability  at  the  magnetopause  (Dungey  [1955],  Parker  [1958],  and 
Southwood  [1968])  may  be  observed  on  the  ground.  However,  those 
waves  should  be  largely  modified  in  their  modes  in  transimission 
across  the  electric  current  layer  of  the  ionosphere. 

In  the  present  work,  the  disturbance  localization  is  assumed  to 
be  in  the  x-direction  only.  As  a  result,  the  major  magnetic  component 
at  the  ground  surface  is  always  in  the  direction  of  localization  at 
all  the  ULF  frequencies  analyzed,  so  that  linear  polarization  occurs 
on  the  ground.  Thus,  the  wave  polarization  of  the  magnetic  field  on 
the  ground  is  quite  sensitive  to  the  localization  of  the  disturbance. 

Consequently,  the  polarization  of  the  wave  disturbances  generated 
in  the  magnetosphere  is  difficult  to  infer  from  the  magnetic  polarization 
observed  on  the  ground  unless  the  localization  of  the  disturbance  is 
precisely  known.  For  example,  when  a  disturbance  is  localized  in 

-Ax 

two  directions  and  is  of  the  form  6  •  £,  °  an  elliptically 

polarized  wave  may  possibly  result  on  the  ground.  Dungey  and 
Southwood  [1970]  have  stressed  the  importance  of  the  neutral 
atmosphere  and  of  the  localization  of  the  disturbance  in  determining 
the  modification  of  hydromagnetic  waves  transmitted  through  the 
ionosphere.  The  modification  of  the  wave  character  by  different 
localizations  and  incident  wave  polarizations  must  therefore  be 
carefully  examined. 

rrom  the  present  analysis,  a  vertical  magnetic  component  at 
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the  Earth's  surface  is  found  to  be  associated  with  all  localized  disturbance 


cases.  Geomagnetic  micropulsations  observed  in  the  auroral  zones 
always  contain  a  vertical  magnetic  component.  The  inference  may  then 
be  made  that  the  geomagnetic  micropulsations  observed  at  the  auroral 
zone  stations  are  originally  localized  disturbances. 

As  Heacock  [1970]  and  Tartaglia  [1970]  have  shown,  the  wave 
polarization  observed  in  the  auroral  zone  depends  on  the  location  of 
the  observatory  relative  to  the  center  of  the  disturbance.  This  fact 
also  suggests  that  the  disturbances  occurring  in  the  auroral  zone  are 
localized.  Thus,  to  discuss  the  observed  wave  polarization  in  terms 
of  plane  waves  in  a  uniform  plasma  is  not  so  meaningful. 

4.36  Comparision  of  Localized  Disturbance  and  Plane  Wave  Cases 
at  1/300  c/s 

In  contrast  with  the  localized  disturbance  the  plane  wave  case  is 
isotropic  above  250  km.  Both  the  slow  and  fast  wave  quartic  roots  have 
the  same  values  above  2  50  km.  as  the  slow  wave  quartic  roots  of  the 
localized  disturbance.  Below  this  level  the  plane  wave  case  is  also 
anisotropic  but  the  behavior  of  the  plane  wave  quartic  roots  is  different 
from  the  localized  disturbance  quartic  roots  in  the  lower  ionosphere 
(See  Figure  4.23).  The  medium  character  is  summarized  in  Table  4.4. 

The  quartic  roots  indicate  the  slow  mode  for  the  plane  wave  case 
is  propagating  throughout  the  ionosphere  whereas  the  fast  mode  is 
propagating  down  to  about  250  km . ,  below  which  the  fast  mode  is 
evanescent.  In  contrast,  for  the  localized  disturbance  case  the  slow 
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TABLE  4.4  Medium  characteristics  at  1/300  c/s  for  a  plane  wave. 


mode  is  propagating  down  to  about  80  km.  where  it  becomes  evanescent 
and  the  fast  mode  is  evanescent  throughout  the  ionospi  3re. 

From  the  dispersion  relation  of  section  2 .2  the  quartic  roots  for 
the  plane  wave  case  above  800  km.  are  given  approximately  by 

for  both  the  slow  and  fast  waves.  Because  the  disturbance  factor  /  is 
zero,  both  modes  <  e  propagating  .  Between  250  km.  am  450  1  i.  the 
slow  and  fast  wave  quartic  roots  are  both  given  approximately  by 

The  effect  of  charge  exchange  on  the  Booker  quartic  roots  is  shown 
in  Figure  4.24  .  Charge  exchange  of  0+  with  O  significantly  enhances 
the  real  part  of  the  Pedersen  conductivity  in  the  region  200-500  km.  and 
causes  this  region  to  be  metallic,  ie. ,  the  real  and  imaginary  parts  of 
the  quartic  roots  are  equal.  If  charge  exchange  is  neglected  this  metallic 
region  occurs  between  200  km.  and  2  50  km. 

Although  charge  exchange  considerably  enlarges  the  imaginary 
parts  of  the  quartic  roots,  and  hence  absorption,  in  this  region  the 
wavelength  is  so  large  that  the  ionosphere  is  still  almost  completely 
transparent  to  the  disturbance.  Also  at  this  frequency  the  coupling 
coefficients  are  larger  than  the  quartic  roots  and  are  the  dominant 
factor  in  determining  the  final  form  of  the  wave  fields. 
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Booker  Quartic  Roots 

Effect  of  charge  exchange  on  the  Booker  Quartic  Roots 


The  neutral  atmosphere  behaves  as  free  space  to  the  plane  wave 
and  as  an  elastic  evanescent  region  to  the  localized  disturbance.  The 
interior  of  the  Earth  is  isotropic  for  the  plane  wave  and  the  localized 
disturbance,  both  of  which  have  the  same  values  for  the  quartic  roots. 

For  the  plane  wave  the  coupling  coefficients  are  the  same  for 
both  the  slow  and  fast  modes.  Also,  the  plane  wave  coupling 
coefficients  are  the  same  as  the  coupling  coefficients  of  the  localized 
slow  mode  except  that  the  mode  coupling  and  mode  reflection 
type  coupling  terms  are  absent  from  the  plane  wave  case  (See  Figure  4.2r>). 
Consequently  reflection  type  coupling  is  dominant  for  the  plane  wave 
case;  mode  coupling  and  mode  reflection  type  coupling  are  a  result 
of  the  localized  nature  of  the  disturbance. 

The  mixing  ratios  reveal  that  the  slow  mode  plane  wave  is  transverse 
magnetic  at  4000  km.  and  inside  the  Earth.  This  fact  is  also  borne  out 
by  the  wave  fields.  The  localized  slow  mode  was  seen  to  change  from 
transverse  magnetic  at  4000  km.  to  transverse  ef^ctric  inside  the  Earth. 

The  incident  fast  mode  is  transverse  electric  at  all  altitudes  for  the 
plane  wave  and  localized  disturbance  cases. 

The  wave  fields  for  the  plane  wave  (Figures  4.2*'  and  1.2'/)  reveal 
that  the  z  components  of  the  electric  and  magnetic  fields  are  zero. 
Therefore  the  plane  wave  will  not  induce  any  surface  charge  at  the  h  >**  >m 
of  the  lower  ionosphere  and  at  ground.  Also,  the  electric  field  of  the 
slow  mode  plane  wave  is  attenuated  by  only  one  order  of  Magnitude 
between  4000  km.  and  ground  surface  and  the  magnetic  field  is  completely 
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Figure  4.26 


Figure  4.27 


unattenuated  down  to  the  Earth's  surface.  Thus  the  Earth's  ionosphere 
and  atmosphere  are  nearly  transparent  to  the  incident  slow  mode  plane  wave. 

The  electric  field  of  the  incident  fast  mode  plane  wave  is  attenuated 
by  two  orders  of  magnitude  and  the  magnetic  field  is  completely 
unattenuated  between  4000  km.  and  ground  surface.  For  the  localized 
case  the  magnetic  and  electric  wave  fields  of  the  incident  fast  mode 
were  attenuated  by  four  and  five  orders  of  magnitude,  respectively, 
between  4000  km.  and  the  Earth's  surface.  Consequently,  introduction 
of  the  localizing  effect  causes  the  fast  mode  to  be  much  more  severely 
attenuated . 

The  wave  polarization  is  almost  linear  at  every  height  for  both 
the  slow  and  fast  plane  wave  cases.  However,  at  higher  altitudes 
for  the  fast  wave  and  at  lower  altitudes  for  the  slow  wave,  a  slight 
ellipticity  is  discernible.  The  rotation  sense  is  right  handed  for 
the  slow  wave  and  left  handed  for  the  fast  wave.  The  hydromagnetic 
wave  polarization  in  a  collisionless  cold  plasma  should  be  basically 
linear  (See  Stix  [1  962]).  However,  elliptical  polarization  may  result 
when  collisions  are  taken  into  account,  because  of  the  greater  effect 
of  collisions  on  the  dispersion  relation  at  hydromagnetic  frequencies 
than  at  higher  frequencies.  Thus,  the  small  ellipticity  in  the  wave 
polarization  results  from  the  effect  of  collisions. 
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5.0  OVERALL  COMPARISON  OF  THE  DISTURBANCE  PROPAGATION 

CHARACTERISTICS  AT  THE  VARIOUS  FREQUENCIES  INVESTIGATED 

In  Tables  5.1  and  5.2  are  found  some  of  the  physical  parameters 
of  the  incident  slow  and  fast  modes,  respectively,  at  5000  km.  for 
all  the  frequencies  considered  in  this  investigation.  The  propagation 
angle  0  is  measured  with  respect  to  the  magnetic  field  line,  i.e., 
from  the  vertical.  The  quantity  &  is  the  skin  depth  measured  in  km. 
and  is  given  by 


where  q^  is  the  imaginary  part  of  the  Booker  quartic  root  and  Re  is  the 
Earth's  radius  in  km.  The  wavelengths  \  x  and  Az  in  the  x  and  z 
directions,  respectively,  are  given  by 

• 

\  zrrRt 

V  ~T~ 

\  .  lUz 

1  “Vr 

where  is  the  Fourier  number  and  qr  is  the  real  part  of  the  Booker 
quartic  root.  Notice  that  at  the  lower  frequencies  is  vety  large, 
i.e.,  about  1.6  Re  at  1/3  cps.  and  160  Re  at  1/300  cps .  for  the 
incident  slow  wave. 

Reflection  coefficients  for  the  various  incident  slow  and  fast 
mode  disturbances  are  tabulated  in  Table  5.3.  The  coefficients  are 
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TABLE  5.1  Disturbance  parameters  for  slow  moHe  at  5000  Km. 


TABU,  5.2  Disturbance  parameters  for  fast  mode  at  5000  Km. 


defined  in  terms  of  the  mixing  ratios  at  5000  km.  For  example,  R  ^ 
means  the  ratio  of  the  upgoing  fast  wave  to  the  downgoing  slow  wave  at 

5000  km. ,  where  r  is  tin  modulus  of  the  ratio  and  ©  is  the  phase  angle. 

i  6 

Thus  R  =  re 

At  100  cps.  the  incident  slow  mode  is  almost  entirely  reflected 
due  to  the  ion  resonance  which  is  encountered  at  1070  km.  At  3.0 
cps.  the  incident  fast  wave  is  well  reflected  whereas  at  1/3  cps.  the 
incident  slow  wave  is  well  reflected.  At  1/300  cps.  the  incident 
slow  and  fast  plane  waves  are  almost  completely  reflected;  Tables 
5.4  and  5.5  reveal  that  the  transmission  coefficients  for  these  waves 
are  nearly  unity  at  ground  level  so  that  most  of  the  reflection  takes 
place  at  the  Earth's  surface. 

The  transmission  coefficients  for  the  electric  and  magnetic  wav^. 

fields  of  the  fast  and  slow  waves,  respectively,  arc  fend  in  Tables  5.4 

r 

and  5.5.  As  an  example  the  transmission  coefficient  TxZ  is  defined 
as  the  ratio  of  the  magnitude  of  the  electric  field  z  component  at  ground 
level  to  the  magnitude  ot  the  electric  field  x  component  at  4000  km. 

Tjj  is  defined  as  the  ratio  of  the  magnitude  of  the  total  electric  field 
at  ground  level  to  the  magnitude  of  the  total  electric  field  at  4000  km. 

The  same  definitions  apply  to  the  m  ic  field  transmission  coefficients. 

In  general  the  incident  fast  wave  transmission  coefficients  arc 
smaller  than  the  incident  slow  wave  transmission  coefficients  except 
at  100  cps.  where  the  slow  wave  transmission  coefficient''  are  zero 
because  the  ion  resonance  condition  is  encountered  at  1070  km.  and 
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E!  CCTRIC  FIELD  MAGNETIC  FIELD 

TRANSMISSION  COEFFICIENTS  TRANSMISSION  COEFFICIENTS 


none  of  the  slow  wave  energy  is  able  to  penetrate  this  barrier.  The 
localized  incident  fast  wave  at  100  cps.  and  the  localized  incident 
slow  wave  at  1/300  cps.  are  very  efficient  in  storing  electric  field 
energy  in  the  neutral  atmosphere.  The  transmission  coefficient  Tj? 
of  the  former  is  3.5  times  larger  than  T^  for  the  latter.  Also,  the 
incident  slow  wave  at  1/3  cps .  is  very  efficient  at  storing  magnetic 
field  energy  below  the  F  region  due  to  energy  being  reflected  back  and 
forth  between  the  Earth  and  the  conducting  wall  appearing  between 
200  and  600  km.  as  a  result  of  charge  exchange  of  0+  with  O,  which 
enhances  the  Pedersen  conductivity. 

In  Table  5.6  a  comparison  is  made  of  magnetic  field  transmission 
coefficients  obtained  by  various  investigators  for  two  different  wave 
frequencies.  For  the  present  work  L  and  R  correspond  to  the  incident 
slow  and  incident  fast  waves,  respectively;  these  are  propagating 
at  16°  and  17°,  respectively,  at  3.0  cps.  and  at  0.0°  for  both  modes 
at  1/300  cps.  The  results  of  Altman  and  Fijalkow  [1968]  were  obtained 
for  vertical  propagation  but  with  the  magnetic  field  inclined  at 20° 
from  the  vertical .  The  results  of  E.O.  Field  [1963]  were  obtained  for 
vertical  polar  propagation;  L  corresponds  to  the  ordinary  mode  and  R 
corresponds  to  the  extraordinary  mode.  Notice  that  the  L  (R)  results 
of  the  present  work  at  3.0  cps.  correspond  more  closely  to  the  R 
(L)  results  of  Field.  This  is  due  to  the  fact  that  the  incident  slow 
(fast)  wave  reverses  polarization  and  propagates  through  the  ionosphere 
primarily  in  the  R  (L)  mode.  Altman  and  Fijalkow  also  found  that  energy 
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Freq 

(c/s) 

INVESTIGATOR 

L 

R 

3.0 

Shaeffer  &  Inoue 

2.25xl0~2 

-3 

1.06x10 

E.C .  Field 

2xl0~3 

8.5xl0"2 

Altman 

and 

Fijalkow 

2.5xl0-2 

1/300 

Shaeffer  &  Inoue 

r 

9.81X10"1 

g.zgxio-1 

Greifinger 

and 

Greifinger 

2x10° 

TABLE  5.6  Comparison  of  magnetic  field  transmission 

coefficients  obtained  by  var;..us  investigators 
for  two  different  wave  frequencies. 
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was  being  fed  from  L  to  R  in  transmission  through  the  ionosphere. 

Greifinger  and  Greifinger  also  considered  polar  propagation  and 
found  somewhat  larner  transmission  coefficients  than  what  the  present 
investigation  reveals . 

Campbell  [1967]  compiled  attenuation  factors  computed  by  various 
investigators  (Akasofu,  1965;  Francis  and  Karplus,  1960;  Greifinger  and 
Greifinger,  1965a;  Karplus  et.  al.,  1962;  Prince  et.  al.,  1964; 

Wentworth,  1964)  and  found  that  Pci  signals  were  reduced  to  about 
0.5  -  0.7  at  0.5  cps .  and  0.2  -  0.5  at  1.0  cps .  of  the  initial  amplitudes 
at  the  top  of  the  ionosphere  in  transmission  through  the  daytime 
ionosphere  at  midlatitudes.  The  present  work  found  transmission 
coefficients  which  were  much  smaller  at  3.0  cps.  for  the  incident 
slow  and  fast  waves  and  at  1/3  cps.  for  the  incident  slow  wave;  however, 
the  transmission  coefficient  for  the  incident  slow  wave  at  1/3  cps.  was 
greater  than  unity . 

The  discrepancies  of  these  results  with  Campbell's  conclusions 
may  be  attributed  to  several  factors. 

1 .  Previous  works  divided  the  medium  into  thick  layers  with  gross 
properties,  whereas  this  work  considers  a  large  number  of  slabs  in 
which  the  ionospheric  parameters  vary  continuously  with  altitude. 

2.  Incident  plane  waves  were  assumed  in  previous  investigations 
whereas  in  this  work  a  localized  disturbance  was  considered. 

3.  Resonant  charge  exchange  collisions  have  apparently  heretofore 
been  neglected.  This  process  has  been  found  to  be  important  in  determining 
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the  final  wave  fields  in  the  Pci  range. 

4.  Campbell's  results  were  based  on  the  measurement  of  one 
component  only.  Table  5.5  shows  that  this  may  lead  to  completely 

erroneous  conclusions.  Here,  at  3.0  cps.  for  the  incident  slow  mode, 

B  —2 

is  approximately  0.1  but  Tg  is  2.25  x  10.  Thus  three  component 

measurements  are  in  general  necessary. 

In  Tables  5.7  and  5.8  are  given  the  magnitudes  of  thn  electric 
and  magnetic  wave  fields,  predicted  from  the  full  wave  calculations, 
at  various  altitudes  for  typical  observed  values,  given  in  the  doubly 
ruled  boxes,  at  some  given  height.  For  the  incident  fast  wave  the 
100  milligamma  magnetic  field  strength  at  the  altitude  hi  (364  km.) 
is  based  on  the  measurements  of  Scarf  et.  al.  [1968].  At  the  other 
frequencies  for  the  incident  fast  wave  typical  ground  observed  magnetic 
field  values  are  chosen,  e.g.,  see  Campbell  [1968],  Also,  for  the 
Pci  and  Pc5  frequencies  the  intermediate  height  hi  is  120  km. 

Notice  that  the  localized  disturbances  at  1/3  and  1/300  cps.  (Table  5.7) 
require  uncommonly  high  magnetic  field  values  at  4000  km.  in  order 
to  observe  a  typical  ground  measurement  of  0.2  and  50  gammas, 
respectively.  Therefore  the  incident  fast  waves  with  the  Fourier  number 
it  equal  to  10 ''TP  at  1/3  and  1/300  cps.  appear  to  be  banished  from 
the  realm  of  reality. 

Consider  the  case  where  the  magnetospheric  plasma  at  4000  km. 
is  moving  perpendicular  to  the  Earth's  magnetic  field  with  a  drift 
velocity  of  0.4  km/sec.  (values  as  high  as  1  km. /sec.  have  been 
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TABLE  5.7  Electric  and  magnetic  wave  field  magnitudes  of  the  incident 
fast  wave  at  various  altitudes.  The  intermediate  height  h. 
is  364  Krn  for  100  c/s  and  120  Km  for  the  other  frequencies. 


observed,  e.g.,  see  Parker  [1968],  p.  21).  Fora  steady  magnetic  field 
-5  2 

value  of  10  wb/m  the  electric  field  generated  at  4000  km.  by 

this  drifting  plasma  is  given  by 

-*> 

E  =  V  x  B 

I  -fc  i  -  3 

or  lEl-  4  x  10  volts/ m . 

If  the  plasma  is  modulated  at  a  frequency  of  1/300  cps.  this  electric 
field  value  is  precisely  what  is  needed  by  either  the  incident  fast  or  slow 
plane  wave  in  order  to  generate  a  50  gamma  PcS  disturbance  at  ground 
level.  In  like  manner  plasma  drift  velocities  between  0.1  and  1 .0 
km. /sec. ,  modulated  at  the  appropriate  frequency,  are  sufficient  for  the 
incident  slow  waves  considered  in  this  investigation  (Table  5.8)  to 
excite  typical  ground  observed  magnetic  field  disturbances  in  the  Pci 
and  Pc5  frequency  ranges. 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

A  localized  disturbance  was  introduced  at  ail  frequencies  con¬ 
sidered  in  this  investigation  and  was  found  to  have  several  effects.  If 
the  localization  was  large  compared  to  the  Pedersen  conductivity,  e.g. , 
at  1/3  and  1/300  c/s,  then  the  incident  fast  wave  was  evanescent  in 
the  upper  regions  of  the  ionosphere.  Also,  the  neutral  atmosphere 
became  an  evanescent  region  and  mode  coupling  became  large  in  the 
ionosphere,  even  for  near-vertical  incidence,  for  both  the  incident 

slow  and  fast  modes.  Localization  resulted  in  a  large  vertical  electric 
field  in  the  neutral  atmosphere  for  the  incident  fast  wave  at  100  c/s  and 
for  the  incident  slow  waves  at  3.0,  1/3,  and  l/30u  c/s. 

Coupling  of  the  various  magnetoionic  characteristic  waves  was 
found  to  be  of  importance  in  determining  final  wave  fields,  in  particular 
at  1/3  c/s  and  1/300  c/s.  In  fact  even  for  vertical  propagation  of  the 
plane  wave  at  1/300  c/s  the  coupling  coefficients  were  found  to  be 
larger  than  the  Booker  quartic  roots  throughout  almost  the  entire  region 
of  interest.  Also,  at  1/300  c/s  the  incident  slow  mode  was  transformed 
from  transverse  magnetic  at  4000  km  to  transverse  electric  in  the  Earth 
as  a  result  of  strong  mode  coupling  which  peaked  around  150  km. 

Inclusion  of  resonant  charge  exchange  collisions  resulted  in  an 
enchanced  Pedersen  conductivity  such  that  the  integrated  Pedersen  con¬ 
ductivity  is  approximately  1.5  times  the  integrated  Hall  conductivity. 
Thus  the  ionospheric  model  of  Jacobs  and  Watanabe  [1962]  for 
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hydromagnetic  wave  propagation  which  neglected  the  Pedersen  region 
must  be  revised  to  include  two  conducting  layers.  Magnetic  energy  is 
stored  in  the  ionosphere  as  a  result  of  reflection  from  the  Pedersen 
layer,  viz  for  the  incident  fast  wave  at  100  c/s,  for  the. incident  slow 
and  fast  waves  at  3.0  c/s,  and  for  the  incident  slow  wave  at  1/3  c/s. 
Also,  at  1/3  c/s  magnetic  wave  cneigy  for  the  incident  fast  wave  is 
stored  between  about  200  km  and  600  km,  where  the  upper  limit  was 
determined  by  the  localizing  factor,^ which  caused  the  incident  fast 
wave  to  be  evanescent  above  600  km. 

Charge  exchange  also  enhances  the  imaginary  part  of  the  quartic 
rents  so  that  some  absortion  of  wave  energy  takes  place  in  the  Pc  1 
frequency  range  between  approximately  150  km  and  250  km.  The  absorb- 
tion  effect  is  not  as  clear  in  the  Pc  5  frequency  range  where  the  coupling 
coefficients  are  everywhere  larger  than  the  quartic  roots. 

On  the  Earth's  surface  the  observable  wave  field  is  a  residual 
between  the  Incoming  and  reflected  waves  and  is  usually  small 
because  of  almost  perfect  reflection.  Also,  this  residual  is  sensitive 
to  the  conductivity  of  the  Earth's  surface.  Therefore,  the  assumption 
of  a  perfectly  conducting  Earth  is  not  appropriate  for  the  purpose  of 
comparision  of  the  wave  characteristics  on  the  ground  and  at  the  top 
of  the  ionosphere.  Reflection  at  the  Earth’s  surface  can  affect  the  , 
wave  fields  in  the  hydromagnetic  frequency  range  throughout  the  entii^p 
medium  being  considered.  The  reflection  strongly  depends  on  the  form 
of  the  disturbance  as  well  as  on  the  conductivity  of  the  ground  surface 
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layer.  Thus,  it  is  worthwhile  to  investigate  the  resultant  wave  field  from 
the  full  wave  solution  by  changing  the  ground  surface  conductivity  and 
the  functional  form  of  the  disturbance.  The  surface  conductivity  data 
obtained  by  the  Dallas  geomagnetic  station  and  used  in  the  present 
analysis  are  perhaps  somewhat  high  compared  to  data  obtained  by 
stations  located  in  areas  of  different  geological  structure. 

The  incident  fast  waves  with  a  disturbance  factor^  value  of  io<\fz7 
at  1/3  and  1/300  c/s  are  deemed  incapable  of  being  observed  on  the 
ground  in  view  of  the  extremely  high  field  values  necessary  at  4000  km 
in  order  to  generate  typical  ground  observed  values. 

At  4000  km  altitude  plasma  drift  velocities  between  0.1  and  1 .0 
km/ sec,  modulated  at  the  appropriate  frequency,  are  sufficient  for  the 
incident  slow  waves  considered  in  this  investigation  fTable  5.8)  to 
excite  typical  ground  observed  magnetic  field  disturbances  in  the  Pc  1 
and  Pc  5  frequency  ranges . 

At  3.0  c/s  (Pc  1)  the  incident  slow  (fast)  mode  was  found  to 
undergo  polarization  reversal  at  1000  km  and  to  propagate  through  most 
of  the  ionosphere  polarized  in  the  right  deft)  hand  sense,  although  the 
polarization  at  ground  level  was  linear.  In  view  of  the  fact  that  the 
incident  fast  wave  is  highly  attenuated,  this  conclusion  seems  to  be 
consistent  with  Campbell's  (1967]  hypothesis  that  ground  observed  signals 
have  characteristic  R  wave  absorption  in  the  ionosphere. 
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6.2  Recommendations 


In  order  to  determine  more  fully  the  nature  of  electromagnetic 
disturbances  in  the  ULF  and  ELF  frequency  range  the  following  investi¬ 
gations  are  considered  to  be  necessary. 

1 .  Narrow  band  measurements  of  electric  and  magnetic  wave  fields 
in  the  ELF  frequency  range  at  satellite  altitudes. 

2.  Polarization  and  wave  field  magnitude  measurements  of  micro- 
nulsation  disturbances  at  satellite  altitudes. 

3.  Further  full  wave  analysis  at  different  frequencies  and  varying 
ionospheric  conditions. 

4.  Simultaneous  observation  of  satellite  and  ground  based  signals. 
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1.0  INTRODUCTION 


1 . 1  Purpose 

A  study  of  the  electrically  conducting  properties  of  the  ionosphere 
is  of  particular  interest  in  the  problem  of  wave  propagation  through  the 
ionosphere  because  the  dispersion  relation  may  be  written  in  terms  of 
the  various  elements  of  the  conductivity  tensor.  This  is  effected  by 
writing  Stix's  quantities  P,R,L,S,  and  D  in  terms  of  the  ionosphere's 
electrical  conductivity.  As  a  result,  the  height  dependent  characteristics 
of  the  Booker  quartic  roots  may  be  examined  in  terms  of  the  height 
dependent  behavior  of  the  Hall,  Pedersen,  and  parallel  conductivities. 

1 . 2  Historical  Review 

The  electrical  conductivity  of  the  ionosphere  has  been  of  interest 
to  geophysicists  ever  since  Balfour  Stewart  [1882]  suggested  that  the 
ionospher"  must  contain  an  electrically  conducting  region  in  order  to 
explain  geomagnetic  quiet  daily  variations .  Stewart  theorized  that  the 
solar  heating  of  the  upper  atmosphere  produced  convective  motions  of 
conducting  air  across  the  Earth's  magnetic  field  lines  generatinq  electrical 
currents  by  mearu  of  a  dynamo  action. 

Stewart's  theory  was  placed  on  firmer  ground  by  Schuster  [1889, 
1908]  who  used  a  spherical  harmonic  analysis  to  show  that  the  quiet  daily 
variation  could  originate  above  the  earth  and  also  gave  a  mathematical 
formulation  to  Stewart's  dynamo  theory.  Kenneiiy  [1902]  and  Heaviside 
[1902]  independently  inferred  the  existence  of  an  electrically  conducting 
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where  the  k  and  j  subscripts  re£er  to  positively  and  negatively  charged 
particles ,  and  vice  versa .  The  quantity  e  is  the  electronic  charge  and 
is  considered  positive.  The  quantity is  positive  or  negative  depending 
on  whether  the  charge  is  positive  or  negative.  The  electric  and  magnetic 
fields  are  represented  by  E  and  B,  respectively.  The  part'  :le  velocity, 
density,  mass,  and  mass  density  are  given  by  "v ,  N,  m,  and  ^  , 
respectively.  *V  ^  refers  to  the  collision  frequency  for  positively  and 
negatively  charged  particles.  The  subscripts  e,  i,  and  r.  refer  to  electrons, 
ions,  and  neutral  particles. 

The  terms  involving  temperature  have  been  neglected.  Thus  the  cold 
plasma  approximation  is  made  and  only  two  modes  will  be  obtained  from  the 
disperson  relation.  The  magnetosonic  mode  appears  only  when  temperature 
effects  are  included.  The  gravitational  force  is  neglected  and  has  no 
effect  in  the  micropulsation  range.  (For  a  further  discussion  of  this 
latter  point  sec,  for  example,  Hultqvtst  [1966H  The  ionosphere  is  con¬ 
sidered  to  be  a  two  component  plasma  so  that  Ne  =  Nj. 

Also,  the  neutral  particles  may  be  considered  to  be  stationary  av 
least  as  long  as  the  following  condition  is  fulfilled: 

V  »C0 

where 
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and  U)  is  the  angular  wave  frequency  of  the  disturbance.  A  plot  of 
the  quantity  *V  vs.  altitude  is  shown  in  Figure  1.1.  From  this  figure  V 
is  seen  to  be  significantly  large  only  around  the  F  layer  maximum  and  even 
in  this  region  "V  is  sufficiently  small  for  the  neutral  particle  motion 

_3 

to  be  neglected  for  angular  frequencies  greater  than  2.3  x  10  rad/sec. 

1 .4  Definition  of  the  Average  Ion 

In  general  the  equations  of  motion  are  coupled  through  collisions . 
However,  coupling  of  the  equations  of  motion  is  found  (see  section  4.2  of 
chapter  4)  to  be  unimportant  in  determining  the  transmission  properties  of 
a  disturbance  originating  at  5000  km.  and  propagating  to  the  Earth. 
Therefore  the  equatic  s  of  motion  are  assumed  to  be  uncoupled,  i.e.,  the 
neutral  particle  velocity  is  set  to  zero  and  the  term  involving  the  difference 
in  the  electron  and  ion  velocities  is  replaced  by  the  velocity  of  the 
particle  whose  equation  of  motion  is  being  considered. 

Then  the  equation  of  motion  for  the  jth  ion  species  is  given  by 


Y"j£U6Vf+V®)-TiVa  <1"4' 

where  V  j  is  the  appropriate  collision  frequency  for  the  jth  positive  ion 
with  all  other  particles,  which  we  see  from  the  original  equation  of  motion 
to  be  given  by 
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where  'V  is  the  collision  frequency  of  the  jth  positive  ion  with  the  nth 
T  jn 

neutral  particle,  due  to  either  resonance  charge  exchange  or  polarization. 
Linearizing  and  summing  all  the  transport  equations  for  all  the  ion  species: 


To  simplify  matters,  we  want  to  reduce  to  one  equation  of  the  form 


it.)  - 


The  assumption  is  made  that 


V 

Vyi 

i.e. ,  all  the  ion  species  are  characterized  by  the  same  macroscopic  drift 
velocity.  Then  we  obtain  the  result  that 


and  the  collision  frequency  is  given  by 


v.- 


ly"3(^^vd"V] 
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where 


and 


N.M; 


V  = 

w\ 


The  quantities  Nj  and  represent  the  "average  ion"  mass  and  total  ion 
particle  density.  Is  the  total  averaged  ion  collision  frequency 

including  collisions  with  electrons  and  neutral  particles. 

In  the  same  fashion  the  linearized  electron  equation  of  motion 
is  given  by 


-  t  ~  -ehifc(  £  +  VeX  -  ^wvtVeVt 

where 

,w) 

-  ^  +  Vg.c 

i.e.,  “Vq  is  the  >.otal  electron  collision  frequency,  'Ven^  is  the  collision 
frequency  for  electrons  with  the  kth  neutral  species,  and  *yej  .  .he  collision 
frequency  of  electrons  with  the  jth  ion  species. 

We  have  assumed  throughout  this  discussion  that  all  the  positive 
ions  are  singly  charged. 
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2.0  COLLISIONS 


2.1  Types  of  Collisions 

Table  2.1  shows  the  types  of  collisions  which  are  included  in  the 
collisional  terms  of  the  equations  of  motion.  Resonance  charge  exchange 
collisions,  heretofore  apparently  not  included  in  the  equations  for  wave 
propagation  through  the  ionosphere,  are  found  to  be  an  important  process 
and  expressions  given  by  Banks  [1966b]  are  utilized.  For  electron-ion 
and  ion-electron  collisions  the  Fokker- Planck  equation  is  used  to  give  a 
slightly  improved  collision  frequency  (see  section  2.4)  over  previous 
expressions , 


Table  2.1 


'^\^olliding 

IncideM^C^ter 

Particle 

Neutral  Particle 

Electron 

Pos .  Ion 

Charge 

Exchange 

Polarization 

Empirical 

Fokker 

Planck 

Fokker 

Planck 

Electron 

X 

X 

Positive  Ion 

X 

X 

X 

In  general,  electron-electron  and  ion-ion  collisions  may  be  ignored  (see 
discussion  in  section  2.4). 

2.2  Electron-Neutral  Collisions 

The  collision  frequency  in  the  friction  term  of  the  equations  of 
motion  has  long  been  a  matter  of  controversy.  In  magnetoionic  theory 
[Appleton,  1932;  Ratcliffe,  1959]  the  average  number  of  collisions  which 
an  electron  makes  per  unit  time  is  assumed  independent  of  electron  velocity. 
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Experimental  evidence  indicates  that  the  collision  frequency  is  velocity 
dependent  and  that  this  dependence  varies  with  the  type  of  plasma.  Phelps 
and  Pack  [1959]  have  shown  that  the  electron  collision  frequency  in 
nitrogen  gas  is  proportional  to  the  square  of  the  velocity.  This  experi¬ 
mental  evidence  motivated  Sen  and  Wyller  [1960]  to  re- formulate  the 
Appleton- Hartree  magnetoionic  formulae  taking  into  account  the  velocity 
dependence  of  the  electron- neutral  collision  frequency.  They  assumed  the 
neutral  particle  velocity  distribution  remained  Maxwellian  under  the 
influence  of  an  alternating  electric  field  and  a  steady  magnetic  field  and 
then  expanded  the  electron  velocity  distribution  (using  the  Chapman- 
Enskog  method)  to  first  order,  retaining  the  zero  order  isotropic  part 
and  the  first  order  anisotropic  part.  The  collision  integral  on  the  right 
hand  side  of  the  Boltzman  equation  was  assumed  to  have  the  form  given 
by  Chapman  and  Cowling  [1958,  p.  350]  which  contains  the  collision 
frequency.  Sen  and  Wyller  then  proceeded  to  develop  the  conductivity 
tensor  and  indices  of  refraction  in  terms  of  C  script  integrals  which  have 
been  tabulated  by  Dingle,  Arndt,  and  Roy  [1957],  These  integrals 
allow  one  to  take  into  account  any  functional  form  of  the  collision 
frequency  which  can  be  expressed  as  a  power  series  in  the  electron 
velocity.  Sen  and  Wyller  found  that  their  absorption  factors  could  differ 
by  as  much  as  100  per  cent  from  the  results  of  the  Appieton-Hartree 
formulae. 

The  electron-neutral  collision  frequency  expressions  utilized  in 
the  present  work  were  adopted  from  the  results  of  Ranks  [1966a] .  Banks 
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made  a  study  of  elastic  collisions  between  electron  and  neutral  gases 
characterized  by  separate  Maxwellian  velocity  distributions  and  then 
defined  collision  frequencies  for  electrons  colliding  with  N2,  C^,  O,  H, 
and  He  in  terms  of  empirical  momentum  transfer  cross  sections,  except  for 
the  cross  sections  of  O  and  H  which  were  derived  from  theoretical  cal¬ 
culations  of  scattering  phase  shifts. 

These  collision  frequencies  given  by  Banks  are  as  follows: 

Table  2.2 

Neutral  gas  Collision  frequency  Electron  temperature  range 

n2  V  =  Nn2  (1.44X10-11)  Te  232°K  $  Te  $1160°K 

V  =  Nn2  (1.41X10-11) 

X  [l-(3.13xl0‘3)  Tg]  T^  1160°K  ^  Te  ^  11,600°K 

°2  V  =  No2(1.37x10"10)  [1+(2.86x10"2)T^]T^ 

232°K  <  T  <  l l , 600°K 
e 

O  V  =  N  [(2.12*0.623)  x  10"10]  T^  T_  <  4000°K 

H  V=  Nu  (3 .41xl0~9)  [l-  (9 .  Oxio”5)  T  1  T* 
h  e  e 

1 50°K  ^  Te  $  5000°K 

He  NHe  [(3.4&±0.373)  x  lO-10]^  Te*58,000°K 

-3 

where  the  neutral  particle  densities  are  measured  in  cm  . 

For  collisions  of  electrons  with  argon  the  hard  sphere  model  is 
used  and  the  collision  frequency  is  given  by 

■V  *  Na  ft.37xlO‘ll)T^ 
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In  general  theoretical  methods  do  not  seem  to  yield  accurate  values 
for  momentum  transfer  cross  sections  so  that  Banks  relied  as  much  as 
possible  on  experimentally  obtained  values .  His  collision  frequencies 
have  been  incorporated  into  the  equations  of  motion  by  utilizing  equations 
(1-7).  Perhaps  a  more  accurate  method  would  have  been  to  use  the  velocity 
dependent  momentum  transfer  cross  sections  given  by  Banks  in  the  equations 
developed  by  Sen  and  Wyller  [I960], 

2.3  Ion- Neutral  Collisions 

2.31  Resonant  Charge  Exchange 

Banks  [1966b]  found  that  resonant  charge  exchange  between  an 
ion  and  its  parent  neutral  molecule  can  significantly  enhance  the 
magnitude  of  the  momentum  transfer  cross  section.  An  example  of  this 

l  4* 

type  of  process  is  O  +  O  - ►  O  +  O  .  T  iere  is  a  transfer  of  charge 

between  two  particles  but  each  approximately  retains  its  original 
kinetic  energy.  Therefore,  even  though  the  ion  has  changed  its  identity 
the  reaction  is  basically  elastic.  This  reaction  provides  a  means  of 
transforming  energetic  ions  into  energetic  neutral  particles .  The 
enhancement  of  the  momentum  transfer  cross  section  is  a  consequence  of 
the  conversion  of  a  grazing  collision  into  an  apparently  direct  impact 
resulting  in  significantly  large  back  scattering,  as  illustrated  in  the 
following  diagram. 
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The  dotted  lines  indicate  the  apparent  backscattering. 

The  resulting  charge  exchange  collision  frequency  for  momentum 
transfer  given  by  Banks  is  changed  to  the  form 

^C.E.  =  (5.93x10  )£ji  ( r+T)^  [E  +  D  •  Ln  fTj  +  Tn) ] 

^  1  0 

-3 

where  Nn»  Mn#  T^,  and  T^  are  the  neutral  particle  density  in  cm  , 

o 

neutral  particle  mass  in  kg,  and  neutral  and  ion  temperatures  in  K, 
respectively.  The  logarithm  is  to  the  base  e  and  the  coefficients  D  and 
E  are  given  in  the  table  below. 

Table  2.3 

Molecule  D  E 


H 

0.144  x  10 

-4.14  x  10 

O 

0.0824 

-3.02 

N 

0.0796 

-2.97 

He 

0.128 

-3.32 

0 

0.0936 

-3.04 

2 

N, 

0.118 

-4.07 

Charge  exchange  sianificantly  enhances  the  total  ion  collision  frequency 
above  150  km  (see  Figure  2.1). 

2.4  Fokker-Planck  Collisions 

A  measure  of  the  ratio  of  large  angle  to  small  angle  coulomb 
collisions  is  the  quantity  (see  Schmidt  [1966],  p.  312): 
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I 

I 

i 


Charge  exchange  effect  on  the  ion  collision  frequency  vs  altitude 


with 


where  \  D  is  the  Debye  shielding  length  and  ♦  is  the  impact  parameter 
for  Coulomb  collisions  when 

X  IVUl 

•to  =  1 

*  4  fAu,v 

M  being  the  reduced  mass  and  u  the  relative  velocity  between  charged 
particles  a  and  b.  For  conditions  existing  in  the  ionosphere,  R  «  1. 

Thus  the  overwhelming  majority  of  particle  encounters  in  the  ionospheric 
plasma  leads  to  small  angle  deflections.  A  Fokker- Planck  treatment  of 
charged  particle  collisions  is  therefore  most  appropriate. 

One  mit,ht  try  to  define  the  small  angle  collision  frequency 
as  the  reciprocal  of  Spltzer's  "slowing-down"  time  or  of  Spitzer's 
"deflection*  time  (Spitzer,  1962,  p.  131).  These  reciprocal  relaxation 
times  are  plotted  in  Figure  3.2  and  Figure  2.3.  Some  of  the  iono¬ 
spheric  parameters  used  in  the  calculation  of  the  various  collision 

t 

frequencies  are  shown  in  Figure  2.4  where  40  .  40..  XL  -  and  XL* 

*  e  ^ 

are  the  electron  and  ion  plasma  frequencies  and  the  electron  and  ion 
cyclotron  frequencies,  respectively.  Spitzer's  relaxation  times  are 
seen  to  be  highly  dependent  on  the  choice  of  the  thermal  velocity  of 
the  test  particle. 
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Componsion  of  Spttzer's  collision  frequency  ond  reciprocal  reloxation  times  for  electrons 


700 


700 


(unoaanunv 


FREQUENCY  (sec-1) 

Plasma,  cycloiron, and  collision  frequencies  vs  altitude  for  electrons  and  ions 


A  more  appropriate  approach  to  the  determination  of  the  Fokker- 
Planck  collision  frequency  is  as  follows.  The  assumption  is  made  that 
the  electron  and  ion  gases  are  both  represented  by  velocity  distributions 
which  are  isotropically  distributed  in  velocity  space  and  are  separated 
by  a  macroscopic  relative  velocity  V.  If  the  velocity  distributions  are 
isotropic  then  when  the  second  moment  of  the  Fokker-Planck  equation  is 
taken,  only  the  friction  term  appears  and  the  diffusion  term  vanishes. 

Also  electron-electron  and  ion- ion  collisions  do  net  appear.  The  velocity 
distribution  for  ions  may  be  approximated  by  a  delta  function  unless  the 
ion  temperature  becomes  extremely  high.  Under  these  conditions  the 
frictional  force  per  electron  is  found  to  be 

f  =  -*v?  •veL£x(  \  +4  +  %  ^ . ) 

where 


This  expression  for  the  electron-ion  collison  frequency  is 

greater  than  the  result  given  by  Spitzer  and  Harm  [1953],  which  is  put 
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into  a  more  convenient  form  by  Shkarofsky  et.  al.  [1966],  and  the 
expression  given  by  Banks  [1966a] . 

In  general  the  quantity  X«1 .  The  electron  and  ion  macroscopic 

drift  velocities  are  \i  and  V-  >  respectively,  and  V  =  V#  ~  V- 

e.  v  &  i 

Then  conservation  of  momentum  requires  that 


Under  the  assumption  that 


it  follows  that 


W«-  A)  . 

M;  e-t 


The  various  collision  frequencies  as  a  function  of  altitude  are 
shown  in  Figures  2.1  to  2.4.  The  Fokker- Planck  correction  becomes 
important  above  150  km.  for  electrons  and  above  1000  km.  for  ions. 
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3.0  CONDUCTIVITY  TENSOR 


3.1  Conductivity  Tensor  Derived  Prom  the  Uncoupled  Equations  of  Motion 
The  uncoupled  equations  of  motion  for  electrons  and  ions  under  the 
influence  of  oscillating  electric  and  magnetic  fields  are  given  by 

i -e(.e 

A\l.  a  V 

S  e'U  *ViX0  -^^V; 

where  B  is  the  Earth's  magnetic  field  and  -ty  and  V.  are  the  ion  and 

v  l  C 

electron  collision  frequencies  defined  by  equations  (1-6)  and  (t-7), 
respectively.  The  conductivity  tensor  is  now  derived  under  the  following 
assumptions: 

1 .  The  equations  of  motion  are  linearized  assuming  the  time 
variation  of  all  quantities  is  of  the  form  £  i,U”C  . 

2.  The  zero  temperature  approximation  is  made,  i.e.,  the  thermal 
motions  of  the  electrons  and  ions  are  ignored. 

3.  The  various  ions  in  the  ionosphere  are  all  represented  by 
an  "average"  singly  charged  ion. 

4.  Negative  ions  are  ignored. 

The  equations  of  motion  are  solved  for  the  electron  and  ion 
velocities  and  then  plugged  into  the  relation  for  the  current  density 


This  result  is  then  compared  with  the  relation 

T  *  W'  E 

to  obtain  the  conductivity  tensor  CT  . 

After  some  algebraic  manipulation  one  obtains  the  result 


-<rH 


where  ,  and  are  the  Pedersen,  Hall,  and  parallel  con¬ 

ductivities,  respectively,  and  are  given  by 


v  iM  +  Jk^L" 

L  \  l  v;*4.n;  v> a] . 


= 


%  - 


eMef 

-Q_t 

A  L 

eNfft  [ 

‘ 

-B 

P  [ 

L  < 

where 


V 


ls*ji 


The  magnetic  field  was  assumed  parallel  to  the  z  -axis  in  this  derivation. 
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Also,  the  vector  cross  product  of  the  velocity  and  magnetic  field  is  put 
in  the  form 
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In  like  manner, 


O 


(«i  v  =  U.5>*eYti'>i 


It  can  be  easily  shown  that 


Therefore 


Then 


P  --  5,  Q.-2.  a.  f 

%  =  ^"0'e^  +  aa;/ix)4 

^  =  £rte(.v_-vt) 

'  +[*u4- BeJillV 


From  momentum  conservation 


V  -  -H  \> 

»c“  W  ^ 

V 

and  therefore 


ai‘y*  ■  2«cy-  =  -  + — 1 

1,1  wx  L  m€  V*.  J 
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Then 


f  =  6ie.P"la;6e-Ye^ 

Now 

y-v/ec 

--  -^(|Hy+^y,2e; 

-YeO+iU-iy^f 

=  -ii[n(i  +  't2eMWtO+'l2;J 

+  ltY^«,-Yei^] 

1 — 5 - 1 

=  -^[Yi''M-Y.V;^Uy;€eJl 

Therefore 

T  =  ^r''[-UY;Vt„-ytV;j 
=  -1-^  n"C'(iVM-YeV;,')e 

wo 

and  the  conductivity  tensor  from  the  ~~up{ed  equations  of  motions  is 

<r«-i^gr-  P  (X^rt-yeV4n) 
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The  dielectric  tensor  is  then 


K-X 


+- 1 


where 


-  V 


The  above  expression  for  the  conductivity  tensor  is  not  very  amenable  to 
solving  explicitly  and  is  therefore  solved  implicitly  with  the  computer. 

3.3  Conductivity  Profiles  for  a  Model  Ionosphere 

A  model  ionosphere  for  midlatitude  and  midday  at  sunspot  minimum 
is  used  to  calculate  the  ionospheric  conductivity  profiles  presented  in  the 
following  sections . 

3.31  Parallel  Conductivity 

Ttgure  3.1  show.?  that  the  Spitzer  arJ  Fokker- Planck  formulas  give 
differs  '  results  I  :r  t'ru  real  part  ct  the  perailei  conductivity  above 

i  jU  kr, .  The  PoKkcr- Planck  case  iu  reduced  b\  a  factor  of  ^TT*  at  ^e 
higher  altitudes.  The  dashed  curves  show  that  the  parallel  conductivity 
is  only  slightly  affected  in  the  range  300-600  km.  by  theincluslon  of 
charge  exchange  effects. 

Figures  3.2  and  3.3  show  the  frequency  variation  of  the  parallel 

conductivity  when  the  Fokker- Planck  and  Spitzer  formulas  are  used, 

respectively.  Above  $00  km.  (£ris  only  slightly  affected  by  varying  the 

-4  o 

frequency  between  10  and  10  c  s.  The  real  part  of  the  parallel  con¬ 
ductivity  is  sianificantly  decreased  above  300  km.  by  increasing  the 

o  1 

frequency  above  !0  cps.  At  700  km.,  for  a  frequency  of  10  cps.. 
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PARALLEL  CONDUCTIVITY  (mhos/meter) 


exchange  and  Fokker-Planck  collisions)  vs  altitude 

Figure  3.2 


GJir  is  approximately  one  order  of  magnitude  smaller  than  0#^- 

o 

for  frequencies  less  than  10  cps. 

3.32  Hall  Conductivity 

Figure  3.4  shows  that  the  real  part  of  the  Hall  conductivity  <T 

Hr 

is  siynificantly  affected  above  120  km.  if  the  charge  exchange  processes 

are  included.  The  Spitzer  and  Fokker- Planck  formulas  give  no  significant 

differences  in  the  magnitude  of  the  Hail  conductivity.  The  frequency 

dependence  of  is  depicted  in  Figure  3.5.  For  frequencies  between 

-4  1 

10  and  10  cps.,  the  real  part  of  the  Hall  conductivity  below  150  km.  is 

-4 

unaffected  by  frequency  variations.  For  frequencies  between  10  and 
o 

10  cps.  is  somewhat  diminished  for  increasing  frequency  above  150  km. 

Also,  actually  becomes  negative  for  the  frequencies  10°  and  101  cps. 
at  the  altitudes  300  and  150  km.,  respectively.  This  transition  is  the  ion 
cyclotron  resonance  effect  In  the  presence  of  collisions  and  is  given  by 
the  condition 


■2- 


with  Y«U) 


i 


where  Y  and  £.<  are  the  usual  Appleton- Hartree  quantities  . 

The  significance  of  the  negative  sign  is  that  the  ions  move  faster 
than  the  electrons.  This  is  due  to  the  fact  that  wave  energy  is  being 
fed  into  the  ion  motion  with  the  result  that  the  ion  gyroradius  incr-  ->s 
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and  enables  the  ions  to  outrun  the  electrons.  Inspection  of  the  imaginary 


part  of  the  Hall  conductivity  shows  that  it  is  larger  than  the  real  part, 
indicating  that  wave  energy  is  being  absorbed.  The  effect  of  ion  collisions 
is  to  reduce  the  efficiency  with  which  the  ion  motion  increases. 

3.33  Pedersen  Conductivity 

Figure  3.6  shows  that  the  real  part  of  the  Pedersen  conductivity 
CTLr  is  significantly  affected  above  130  km.  by  the  inclusion  of  charge 
exchange  effects.  Above  200  km.  CJj^  is  enhanced  one  order  of  magnitude 
or  more.  This  enhancement  is  due  primary  to  charge  exchange  of  O 
with  O.  Although  not  shown  here  the  same  effect  is  also  true  for  the 
imaginary  part  of  the  Pedersen  conductivity  Q ^  and  this  is  very  important 
in  altering  the  index  of  refraction  for  a  disturbance  propagating  through 
the  ionosphere  at  ULF  frequencies  (see  section  4.2  of  Part  II  and  sections 
4.11  and  4.31  of  Part  I)  . 

Figure  3.7  shows  that  0^  is  unaffected  by  frequency  variations 
for  frequencies  less  than  10°  cps .  At  a  frequency  of  101  cps.  C£r  is 
enhanced  slightly  at  all  altitudes  above  130  km. 

3.34  Comparison  of  the  Various  Conductivities 

The  real  parts  of  all  the  conductivity  elements  are  shown  together 
for  comparison  in  Figure  3.8.  The  Hall  conductivity  is  larger  than  the 
Pedersen  conductivity  up  to  about  130  km.  The  maximum  in  occurs 
around  120  km.  A  primary  maximum  in  O'  occurs  at  130  km.  due  to  charge 

±r 

*4* 

exchange  of  O  with  0„  .  A  secondary  maximum  in  0“  occurs  around  200 
2  4 

4* 

km.  due  to  charge  exchange  of  O  with  O.  The  parallel  conductivity  is 
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500 


much  larger  than  the  Hall  and  the  Pedersen  conductivities  throughout  the 
entire  ionosphere. 

Figure  3.9  shows  the  dimensionless  Hall  i5lC  and  Perterc^n  Qlr 

<oeD 

conductivities  at  the  frequency  0.00333  cps.  plotted  on  a  linear  scale  vs 
altitude.  The  integrated  Pedersen  conductivity  (i.e.,  the  area  under  the 
curve)  for  the  F  region  is  1.5  times  the  integrated  Hall  conductivity  for 
the  E  region.  This  represents  a  much  larger  F  region  conductivity  than 
gi.en  by  previous  investigators. 
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4.0  STIX’S  QUANTITIES  IN  TERMS  OF  THE  CONDUCTIVITY  TENSOR 
ELEMENTS 


4.1  Derivation  of  P,R,L,D,  and  S  in  the  Presence  of  Collisions 

If  the  uncoupled  equations  of  motion  are  considered  then  the 
quantity  representing  the  mass  of  the  kth  particle  in  the  absence  of 
collisions  may  be  replaced  by  (1+iZj,,)  for  the  case  when  collisions 
are  included  by  means  of  the  friction  term.  It  is  then  straightforward 


to  show  that  P.R.L.S,  and  D  are  given  by 

R-  .  -  f  vTf* 

L  *  1  -  5T  _J<X— 

*  V* 

S  ■  1/2  (R  +  L)  Dm  1/2  (R  -  I.) 


where 


uv 


y  = 

V  IQ 


*  w 


V  l  +  i** 


"ff  *  plasma  frequency  of  kth  particle 
±1^  *  gyro  frequency  of  kth  particle 


After  some  algebraic  manipulation  these  quantities  may  be  written 


in  terms  of  the  conductivity  tensor  elements  in  the  following  manner 

OU  •  <TT 
R  =  1  +  — JL.  +•  l  — 

C06,  U€, 

i.  -  i  -  .Sl,  +  i  Pi_ 
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P=1+lwt 

s  =  i  +  i, 

W60 

d  = 

wec 

Separating  the  real  and  imaginary  parts 

”  °Hr  *  L  ^rt't 


°L  -  &±x  -V-  ^ 


°n  ~  °//r  v  1  ^ft 


’■U'S&F 


6xr  +•  0\\i 

w  en 


crHr +-(51^ 


‘M 


CUr  - 


ee  J  L  *O60  J 


s  =  j 


£LL  4.  't  j?Ir 
W£«  W£a 


D  =  +  1 


In  the  absence  of  collisions 


-  «Ir  -  °iy  =  © 
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then 


L  =  1  -  ±  C~^ 


s  -  i  -  HjlL 

D  = 

U>£„ 

4.2  The  Effect  of  the  Coupled  Equations  of  Motion  on  Stix's  Quantities 

-3 

In  general,  for  frequencies  above  3  x  10  c/s,  none  of  the 
conductivity  tensor  elements  is  significantly  affected  below  700  km. 
Above  about  1000  km.  the  real  part  of  the  Pedersen  conductivity  and  the 
imaginary  part  of  the  Hall  conductivity  may  be  seriously  affected  by  the 
coupling  of  the  equat’^>ns  of  motion.  However,  the  imaginarv  part  of  th» 
Pedersen  conductivity  and  the  real  part  of  the  Hall  conductivity  are 
more  important  than  the  previously  mentioned  quantities  in  determining 
the  propagation  characteristics  of  a  disturbance  and  are  not  seriously 
affected  by  the  coupling  of  the  equations  of  motion.  Therefore  the 
derivation  of  P,R  ,L,S,  and  D  from  the  uncoupled  equations  of  motions 
will  be  sufficient  for  frequencies  of  interest  in  this  study. 
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PART  III 

A  MODEL  IONOSPHERE 
FOR  MID-DAY  AND  MID- LATITUDE 
DURING  SUNSPOT  MINIMUM 


ABSTRACT 


A  model  ionospnere  for  aaytirne  conditions  at  mid-latitudes  during 
the  period  of  solar  minimum  has  been  constructed  from  existing  experimental 
data.  The  atmospheric  parameters  obtained  from  these  data  include  electron 
and  ion  densities,  electron  and  ion  temperatures,  and  electron  collision 
frequencies  and  cross  sections  for  momentum  transfer. 
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i .  Ion  Densities 


Positive  ion  density  measurements  have  been  extensively  made  up  to 
1000  km.  Few  data  are  available  for  higher  altitudes.  For  negative  ion 
densities,  only  a  few  rocket  measurements  have  been  performed,  but  existing 
data  seem  to  indicate  that  negative  ions  may  be  important  up  to  70  km  during 
the  daytime  and  up  to  90  km  during  the  nighttime  [Sagalyn,  19  65]. 

1  •  1  Ion  Densities  in  the  F  Region 

Positive  ion  densities  in  the  upper  atmosphere  during  solar  minimum  have 
been  reported  Ly  Johnson  [1966],  These  results  are  presented  in  Fig.  1. 

Johnson's  model  is  based  on  two  sets  of  data  which  were  mear.  ared  with 
rud'o  frequency  ion  mass  spectrometers  that  determined  the  relative  ion 
composition.  The  first  data  set  is  for  the  region  90-24C  km  and  was  measured 
from  an  Aerobec  15')  rocket  launched  from  White  Sands,  New  Mexico  at  1634  UT 
(093  I  MSI)  on  February  15,  1963.  The  second  data  set  is  for  the  region 
-lot)- 1200  km  and  was  measured  from  the  Soviet  satellite  Electron  2.  The  latter 
w.is  obta  .,ied  for  the  time  interval  1400  -  1900  hours  LT  and  for  the  latitudes 
10°  -  60°  N  during  the  period  of  February  10-16,  1964. 

Notice  should  be  made  of  the  fact  that  Johnson  derived  the  absolute 
number  densities  for  various  ions  from  the  relative  ion  comp  sition  data  by 
neglecting  negative  ions  and  by  assuming  the  electron  density  was  equal  to  the 
total  ion  density.  To  obtain  this  normalization  Johnson  chose  three  different 
electron  density  Attributions:  ror  90-130  kin  the  results  reported  by  Hr  rdeau 
et.  al.  [l9b6]  obtained  from  a  Nikc-Apachc  rocket,  NASA  14.117,  flown  from 
Wallops  Island,  Va.  at  1207  EST,  November  23,  1964  were  utilized;  the  White 
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Sa-ds  ionosonde  data  reported  by  Holmes  et  al.  [1965]  were  used  for  the  130-240 
km  region;  the  electron  densities  for  the  240-1000  km  region  were  obtained  by 
averaging  the  data  from  a  rocket-borne  dispersive  Doppler  radio  propagation 
experiment  and  data  from  Alouette  satellite,  both  data  sets  being  obtained  500  km 
east  of  Wallops  Island,  Va.  at  0923  EST  {1423  UT)  on  July  2,  1963.  Tohnson 
adopted  electron  densities  for  the  1000-1200  km  altitudes  by  extrapolating  the 
lower  altitude  electron  densities. 

Johnson  presented  no  data  for  the  240-400  km  region.  For  these  altitudes, 
an  interpolation  was  made  between  the  high  and  low  altitude  portions  of  the 
N+  and  0+  curves.  According  to  Johnson,  the  maximum  deviation  required  to 
smoothly  connect  the  three  electron  density  distributions  occurred  at  240  km  and 
was  less  than  25%. 

At  altitudes  above  1200  km  only  a  few  ion  density  measurements  have  been 
made.  Taylor  et,  al.  [19  65]  measured  the  H+  and  He'4  densities  from  1500  km 
to  30,000  km  with  the  positive  ion  spectrometer  on  board  the  OGO-A  satellite 
during  the  period  September  23  through  December  10,  1964.  The  He+  densities 
were  about  1%  of  the  H+  densities  and  thus  were  too  small  to  be  shown  in 
Figure  1.  This  result  implies  that  if  we  intend  to  match  the  data  of  Taylor  et.  al. 
with  Johnson's  model,  then  the  He+  concent"-  on  must  decrease  by  two  orders 
of  magnitude  between  1200  and  1500  km.  Values  for  the  He f  densities  are  given 
in  Table  1  • 

Also  the  total  positive  ion  concentration  N+  was  measured  by  Hanson  [1962] 
with  a  retarding  potential  probe  on  October  4,  1960  during  mid-afternoon.  The 
data  are  given  in  Fig.  1.  Up  to  2400  km  these  values  are  larger  than  those 
given  by  Taylor  et.  al. 
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All  of  the  data  measured  at  solar  minimum  seem  to  indicate  that  H  is  by 
far  the  dominant  ion  above  1200  km  and  that  the  H+  concentration  is  fairly 
constant  between  1200  and  3000  km.  Thus  in  order  to  make  a  consistent  model 
wc  assume  that  H+  is  the  dominant  i  above  1200  km  and  that  the  H+  concen¬ 
tration  exactly  equals  the  electron  density.  We  therefore  extrapolate  Johnson's 
curves  for  the  H+  and  electron  densities  so  that  they  smoothly  connect  at 
about  2000  km  with  the  data  given  by  Taylor  et.al.  This  extrapolation  is 
shown  in  Figure  1  and  values  are  given  in  Table  1 . 

1 . 2  Ion  Densities  in  the  D  and  E  Regions 

The  ionic  composition  ;n  the  lower  ionosphere  is  complicated  by  the  presence 
of  many  minor  constituent  ions.  Fig.  2  shows  these  minor  constituent  ions  to 
be  3  or  4  orders  of  magnitude  smaller  than  the  major  constituent  ions. 

In  addition  to  the  lower  altitude  values  of  Johnson's  model,  two  sets  of 
positive  ion  density  data  are  presented  in  Fig.  2.  The  tirst  set  is  taken  from 
the  results  of  Sagatyn  and  Smiddy,  as  reported  by  Narcissi  [1966],  obtained 
with  an  electrostatic  probe  at  1200  hours  CST,  October  31,  1963  on  the  Nike  Cajun 
AC6.341  launched  from  Elgin,  Florida.  The  second  data  set  was  obtained  by 
Smith  et.  al.  [1967]  with  a  positive  ion  Bennett-type  mass  spectrometer  above 
Wallops  Is.  at  1330  hours  EST  in  May  1962.  Their  measurements  cover  the  mass 
range  7-43  amu . 

Both  sets  of  data  appear  to  be  in  agreement  with  each  other,  except  for  the 
values  of  the  ion  species  24+,  the  values  of  which  given  by  Sagalyn  and  Smiddy 
are  about  3  times  larger  than  those  given  by  Smith  et.  al.  We  choose  the 
Sagalyn  and  Smiddy  data  because  they  were  measured  at  a  time  closer  to  solar 
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minimum.  The  values  of  28+  given  by  Smith  et.  al.  are  adopted  for  the 
aiti  jdes  of  104-112  km.  Then  an  interpolation  is  made  between  these  values 
and  those  measured  by  Sagalyn  and  Smiddy,  whose  results  for  28+  extend  only  up 
to  82  km.  The  interpolated  values  are  indicated  by  the  dashed  line,  We  also 
interpolate  between  Johnson's  values  of  28+  and  those  of  Smith  et.  al. 

The  curves  for  32+  and  30 f  given  by  Smith  et.  al.  are  nearly  parallel  to  the 
corresponding  curves  of  Johnson's  model.  The  32+  curve  of  Johnson  is 
extrapolated  down  to  meet  the  32+  curve  of  Sagalyn  and  Smiddy.  In  the  same 
manner  the  30+  curve  of  Sagalyn  and  Smiddy  is  extrapolated  to  higher  altitudes 
to  join  with  the  30+  curve  of  Smith  et  al .  The  latter  is  then  smoothly  connected 
to  the  30+  curve  of  Johnson's  model. 

In  order  to  provide  a  consistent  model,  we  calculate  the  total  negative  ion 
density  by  subtracting  the  total  electron  density  (to  be  discussed  in  the  next 
section)  from  the  total  positive  ion  density  obtained  by  summing  all  the  densities 
of  the  various  positive  ion  species.  These  results  may  be  compared  with  the 
negative  ion  density  data  reported  by  Sagalyn  [19  65]  and  are  shown  in  Fig.  2. 

Her  measurements  were  made  using  a  two  probe  system  on  a  Nike  Cajun  II  ascent 
at  1200  hours  LT  on  11  March  1964,  Our  negative  ion  density  values  do  not 
agree  with  those  of  Sagalyn  because  our  adopted  altitude  distribution::’  of 
electron  and  positive  ion  densities  are  not  the  same  as  those  reported  by 
Sagalyn.  Taking  our  calculated  values  of  negative  ion  densities  and  plotting 
the  ratio  Ni+/Ne  =  1  +  Mj[/klei  we  obtain  the  curve  shown  in  Fig.  3.  Comparison 
of  our  calculated  ratio  and  the  corresponding  ratio  derived  by  Sagalyn  [1965] 
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from  Nike  Cajun  II  data  shows  that  our  calculated  altitude  distribution  of  Nj/Ne  has 

the  same  shape  above  65  km  as  that  of  Saoalyn,  although  the  negative  ion  density 

values  are  considerably  different.  Below  65  km  Sagalyn's  ratio  of  N./Ne  increases 

Ni/ 

with  decreasing  altitude  whereas  our  calculated  values  of  /  Ne  decrease  to  zero  at 
about  52  km. 

From  Fig.  2  we  sec  that  for  the  region  below  82  km  ions  heavier  than  45  amu 
are  predominant,  with  37  +  ,  19  +  ,  and  30+  being  the  next  most  abundant.  Narcissi 
[1966]  has  identified  these  ions  as  18+  (H20+),  19+  (F^Cr),  28+  or  3i+), 

30 f  (No+) ,  32+'  (0 2)  and  37+  (H5C>2+).  The  heavy  ions  >45  +  ,  as  suggested  by 
Narcisi  [1966],  may  be  higher  hydrates  of  H30+  and/or  other  cluster  10ns. 

Narcisi  argues  that  such  large  quantities  of  hydrates  are  not  expected  from  rocket 
contamination.  Also,  exact  identification  of  individual  negative  ion  constituents 
in  the  lower  ionosphere  has  not  yet  been  made. 

At  an  aititude  of  about.  82  km  a  transition  takes  place  in  the  ionic  composition. 
The  ions  19+  and  37  '  decrease  sharply  while  the  ions  30+  and  32  +  become 
predominant.  At  about  170  km  0+  becomes  the  dominant  ion  species.  As  shown 
in  Fig.  1 ,  a  transition  occurs  at  approximately  1000  km  where  H+  becomes  the 
dominant  ion  and  remains  as  such  at  greater  altitudes.  Values  of  the  adopted 
ton  densities  for  the  D  and  E  regions  are  found  in  Table  2. 
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2.  Electron  Densities 


2 .  1  Electron  Densities  in  the  F  Region 

Several  different  methods  have  been  used  to  measure  electron  densities  in 
the  F  region.  The  results  of  some  of  these  methods  are  presented  in  Fig.  4. 

The  measurements  of  Ulwick  and  Pfister  [19  63]  were  made  with  an  r.f. 
impedance  probe  on  a  rocket  flight  at  0100  hours  EST  on  April  12,  1961.  The 
results  of  Nisbet  et.  al.  [1966]  were  obtained  from  measurements  of  radio  signals 
propagated  between  sections  of  a  Javelin  rocket  flown  from  Wallops  Is.  ,  Va.  on 
May  19,  1965  shortly  after  1600  hours  LT.  The  data  of  Oya  and  Obayashi  [1967  a] 
were  obtained  using  the  sweep  frequency  impedance  probe  technique  on  the 
L-3H-2  rocket  at  1535-1550  hours  JST  on  July  23,  1966.  Jackson  and  Bauer  [1961] 
used  the  Seddon  CW  propa<  tion  technique  on  the  NASA  8.10  rocket  flown  from 
Wallops  Is.  ,  Va.  at  1502  EST  on  April  27 ,  1961 .  Oya  and  Obayashi  [1967  b] 
made  their  measurem  nts  at  1725-1730  JST  on  October  20,  1966. 

These  data  are  rather  spread  but  consistently  indicate  the  presence  of  a 
peak  density  in  the  neighborhood  of  240-320  km.  Also  obvious  from  Fig.  4  is 
the  variability  exhibited  by  the  values  of  the  electron  density  above  about  900  km. 
However,  we  must  remember  that  these  data  were  measured  during  different 
ionospheric  conditions. 

Since  no  negative  ions  have  reportedly  been  measured  at  these  high 
altitudes,  we  are  constrained  to  choose  equal  electron  and  on  density  values 
m  order  to  preserve  charge  neutrality.  We  have  already  adopted  an  ion  density 
distribution  for  the  F  region  (see  Sec  1. 1)  and  must  therefore  select  the  same 
d^iity  distribution  of  electrons.  This  adopted  distribution  is  shown  in  Fig.  4 
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and  comparison  with  the  previously  mentioned  electron  density  data  indicates 
good  agreement. 

2.  2  Electron  Densities  in  the  D  and  E  Regions 

Various  techniques  have  also  been  used  to  measure  electron  densities  in 
the  lower  ionosphere.  The  results  of  some  of  these  methods  are  shown  inFig.  5. 

Smith  Cl 9 65a]  obtained  data  with  Langmuir  probes  during  the  ascent  of  a 
rocket  launched  from  Wallops  Is.  ,  V a.  at  1430  hours  EST  on  February  27,  1963. 

The  data  of  Sagalyn  [1965]  were  measured  with  a  two  probe  system  used  on  the 
Nike  Cajun  II  on  March  ll,  1964  at  1200  hours  LT.  The  data  given  by  Belrose 
et.  al.  [1966]  were  obtained  by  the  partial  reflection  method.  Also,  data 
obtained  from  Nike  Cajun  I  on  October  31  ,  1963  at  1200  hours  LT  with  spherical 
electrostatic  analyzers  are  given  by  Sagalyn  et.  al.  [preprint]. 

To  obtain  electron  density  and  collision  frequency  measurements 
simultaneously,  Salah  and  Bowhill  used  a  particular  radio  propagation  tcchmqur 
by  which  the  transmitting  power  and  polarization  are  controlled  by  the  rocket 
receiver  output  via  telemeter  link.  The  experiment  was  performed  with  a  rocket 
f>red  from  Wallops  Is.  ,  Va.  on  April  12,  1965  at  1214  EST.  Using  the  same 
method  Mechtly  et.  al.  [1967]  obtained  data  from  an  experiment  performed  on 
Nike  Apache  14. 143  launched  from  Wallops  Is.,  Va.  on  April  16,  1964  at  2100  hour 
UT. 

We  choose  curve  1 ,  the  data  obtained  by  Smith,  for  our  model  atmosphere 
because  it  seems  to  be  quite  representative  of  the  various  data  given  in  Fig.  5. 

The  C,D,  and  E  layers  are  clearly  seen  in  Smith’s  electron  density  distribution. 
Ai»ovo  82  km  curve  1  is  in  good  agreement  with  our  adopted  total  positive  ion 
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ELECTRON  DENSITY  (cm"3) 

Figure  5  Electron  density  vs  altitude  in  the  D  and  E  regions. 


density  curve.  Below  about  82  km  Smith's  electron  densities  are  considerably 
less  than  the  positive  ion  densities.  This  is  to  be  expected  because  of  the 
presence  of  negative  ions.  The  negative  ion  density  is  calculated  by 
subtracting  Smith's  electron  density  values  below  82  km  from  the  adopted 
positive  ion  values.  Then  we  extrapolate  the  positive  ion  density  curve  down 
to  52  km  at  which  altitude  the  electron  concentration  is  assumed  to  be  •-'quai 
to  the  positive  ion  concentration.  For  our  model  ionosphere  we  also  assume 
that  the  negative  ions  disappear  above  82  km. 

These  adopted  negative  ion  and  electron  density  values  are  listed  in  Table 
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3 .  Electron  and  Ion  Temperatures 


3 . 1  Electron  and  Ion  Temperatures  in  the  F  Region 

Electron  and  ion  temperatures  in  the  F  region  have  been  measured  by  means 
of  ground  based  radar  backscatter  and  various  types  of  plasma  probes  aboard 
rockets  and  satellites.  A  review  of  these  methods  and  a  comparison  of  the 
results  are  given  by  Evans  [13  63a].  Discrepancies  do  exist  between  the  data 
obtained  by  the  various  methods  and  some  of  these  discrepancies  have  not  yet 
been  resolved.  As  stated  by  Evans,  rocket  and  satellite  measurements  have 
greater  accuracy  and  better  tirru.  and  height  resolution  than  do  radar  measurements. 
However,  up  to  the  present  time,  satellite  data  are  not  sufficient  to  provide  the 
altitude  variations  of  the  temperatures.  Space  vehicles  disturb  the  medium 
surrounding  them.  Extensive  radar  measurements  have  been  made  during  the 
solar  minimum  period  and  over  a  wide  range  of  altitudes.  For  this  reason  we 
choose  the  radar  backscatter  results  for  our  model. 

Ionospheric  backscatter  results  for  solar  minimum  are  given  by  Evans 
[1965a,  1965b,  1967a,  1967b],  These  measurements  were  made  at  the  Millstone 
Hill  radar  observatory  (42.6°N,  71.5°W).  The  results  fur  1963  and  January  1964 
[Evans,  1965b]  are  replotted  in  Fig.  6.  These  data  were  derived  on  the  assumption 
that  only  Q+  ions  are  present  at  all  heights  above  200  km. 

Later  in  1964  Evans  [1967b]  obtained  other  electron  and  inn  temperature 
data  with  different  assumptions.  The  neutral  temperature  above  250  km  was 
determined  from  the  radar  signal  spectra  at  the  altitudes  230,  260,  and  290  km. 

The  neutral  temperature  at  altitudes  below  250  km  was  then  obtained  by  assuming 
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the  exponential  dependence  of  Tn  on  altitude  implied  by  the  1965  CIRA  models 
I  and  Ii  (see  equation  (1).  Evans  [1967b]).  Then  Evans  was  able  to  simultaneously 
determine  Te  3nd  the  ratio  of  (0+)  to  (O^  +  NO+)  ions  for  the  altitudes  135-250  km. 
The  neutral  gas  temperatures  thus  determined  are  as  much  as  200  K°  lower  than 
the  values  given  by  CIRA  1965  for  the  same  solar  flux  (F)  of  10.7  cm  radiation 
averaged  over  5  solar  rotations.  Evan's  neutral  gas  temperatures  disagree  with 
satellite  drag  data  by  9%  [Evans,  19  67b]. 

We  choose  Evan's  data  from  19  63  and  January  1964  for  our  model  for  the 
following  reasons: 

(1)  The  ion  temperature  determined  by  Evans  [1965b]  approximately  equals 

—  _o? 

the  neutral  temperature  given  by  CIRA  19  65  models  I  and  II  (F  ~  65  and  75  x  10 

O 

W/m  c/s,  respectively)  up  to  about  300km.  At  greater  altitudes  Tj  >Tn  This 
result  is  in  fairly  good  agreement  with  theory  [Banks,  19  66’,  Dalgarno  et,  ai .  , 

1967;  Geisler  and  Bowhiil,  1965  a,  1965b]. 

(2)  A  model  for  the  ion  composition  has  already  been  adopted,  as  discussed 

in  sec.  1 ,  and  agrees  with  the  assumption  that  0+  is  the  dominant  ion  above  200  km 

(3)  The  neutral  temperatures  and  corresponding  densities  given  by  the  CIRA 
19  65  model  I  are  adopted  and  there  is  no  need  to  determine  the  neutral  gas 
temperatures  from  radar  signal  spectra. 

Evans  concludes  from  the  data  of  1963  and  January  1964  that  the  altitude 
variations  of  Tfi  do  not  indicate  a  strong  seasonal  dependence.  The  ion 
temperature  shows  even  less  seasonal  variation  than  the  electron  temperature. 

Evans  also  points  out  that  the  temperatures  are  fairly  stable  throughout  the 
daytime.  The  ratio  Tg/Tj  is  greater  than  1.0  throughout  the  F  region  and 
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reaches  a  maximum  value  of  2.0  -  2.6  at  an  altitude  of  about  300  km  soon  after 
dawn  irrespective  of  the  season. 

As  seen  in  Fig.  6,T.  has  a  positive  gradient  at  all  altitudes.  However  there 
remains  some  uncertainty  in  the  gradient  of  Te  above  500  km.  Twelve  different 
profiles  for  Tg  are  depicted.  Nine  of  these  indicate  near-isothermal  conditions 
above  500  km;  three  profiles  have  a  definite  positive  temperature  gradient. 

Theoretical  models  [Banks,  19  66a;  Daigarno  et.  al.  ,  1967;  Geisler  and 
Bowhill,  1965a,  1965b]  predict  that  T0  becomes  isothermal  at  higher  altitudes. 
Evans  [1967a]  has  developed  a  theory  to  explain  the  discrepancy  between  observed 
positive  gradients  and  Geisler  and  Bowhill' s  isothermal  model.  Evans  conduces 
that  this  positive  gradient  and  the  high  protonospheric  temperatures  may  be 
accounted  for  by  the  backward  scattering  of  photoelectrons  from  the  ionosphere. 
This  flux  is  estimated  to  be  about  5  x  10®  eiectrons/cm Vsec  with  mean  initial 
energies  of  14  ev.  However,  the  majority  of  the  cases  presented  in  Fig.  6 
indicates  that  T.is  converging  toward  Te  and  thus  does  not  warrant  invoking 
such  a  theory.  The  twelve  profiles  shown  in  Fig.  6  are  averaged  and  the 
resulting  average  is  shown  in  Fig.  7.  The  values  for  these  average  electron 
and  ion  temperatures  are  given  in  Table  4. 

3 . 2  Electron  and  Ion  Temperatures  in  the  D  and  E  Regions 

Fig  8  shows  the  measurements  of  electron  and  ion  temperatures  made  by 
several  investigators  for  the  D  and  E  regions.  The  data  reported  by  Brace 
et.  al.  [1963]  were  measured  with  a  dumbbell  bipolar  probe  which  was  ejected 
from  the  rocket  vehicle  NASA  6.04  flown  from  Wallops  Is.  ,  Va.  at  1556  hours 
EST  on  March  26,  1961.  Spencer  et,  al.  [J965]  report  results  obtained  with 
a  cylindrical  Langmiur  probe. 
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The  above  data  indicate  a  high  degree  of  variability  of  Tg  in  the  lower 
regions  of  the  ionosphere.  Other  data  [Nagy  et.  al.  ,  1963;  Aono  et,  al.  ,  1962; 
Smith  et.  a.l,  1965b]  not  shown  here  also  substantiate  this  variability  of  Tg. 

Also,  the  comparison  of  ascent  and  descent  data  (see  Fig.  8)  from  individual 
flights  sometimes  reveals  large  differences  in  Te  within  a  time  interval  of  several 
minutes.  This  significantly  large  variation  of  Te  is  not  well  understood  and  makes 
the  task  of  choosing  "average"  data  difficult. 

For  the  sake  of  consistency  we  prefer  to  choose  the  electron  temperatures 
at  these  lower  altitudes  by  extrapolating  the  averaged  Evan's  data  in  Fig.  7 
down  to  120  km  where  we  assume  Te  =  Tn  of  the  CIRA  19  65  model  I  atmosphere. 
This  assumption  is  in  accord  with  theory  [Dalgarno  et.  al.  ,  1967;  Geisler  and 
Bowhill,  1965  a,  1965b].  Values  of  Te  and  Tn  for  this  altitude  range  are  tabulated 
in  Table  4 , 
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4 .  Experimental  Electron  Collision  Frequencies  and  Cross  Sections 
4 . 1  Electron  Collision  Frequencies  for  Momentum  Transfer 
In  Fig.  9  are  reproduced  the  measured  electron  collision  frequencies  for 
momentum  transfer  as  determined  by  various  investigators.  The  values  given  by 
Benson  [1964]  were  determined  with  cross  modulation  experiments  at  College, 

Alaska  on  August  22,  1960,  June  29,  1960,  and  November,  1961.  These  results 
agree,  within  the  probable  error  (+20%),  with  a  formula  Benson  derived  by  adding 
together  the  collision  frequencies  in  N'2  and  O2  and  using  the  pressure  values 
determined  on  a  rocket  flight  in  July  1957.  The  values  of  Jesperson  et.  al.  [1964] 
were  determined  in  August  and  December  of  1962  from  Nike-Cajun  rockets 
(Ferdinand  I  and  II)  which  were  launched  from  Norway  during  auroral  absorption  and 
which  measured  the  relative  amplitude  of  signals  transmitted  from  the  ground  and 
the  phase  relationship  between  signals  of  different  frequencies.  A  marked  difference 
was  obtained  for  the  measurements  made  during  the  two  flights,  viz. ,  the  highest 
values  were  observed  during  summer.  Jesperson  et.  al.  were  not  sure  if  this 
was  a  seasonal  effect  or  only  some  variability  in  the  collision  frequency  values. 

Belrose  and  Burke  [1964]  obtained  collision  frequency  measurements  in 
1961  and  1962  at  Ottawa,  Canada  from  the  recordings  of  the  amplitudes  of  weak 
echoes  of  the  ordinary  and  extraordinary  wave  compenents  reflected  from 
ionization  irregularities.  They  utilized  the  Gardner  and  Pawsey  method  except 
that  the  energy  dependence  of  the  collision  cross  section  was  employed  in  the 
generalized  magnetoionic  formulas.  Belrose  and  Burke  point  out  that  the  observed 
variability  in  collision  frequency  data  may  not  be  due  to  a  seasonal  change  but 
rather  a  variation  with  solar  activity. 
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Hall  and  Fooks  [1965]  determined  the  effective  electron  collision  frequency 
Veii  with  a  radio  propagation  experiment  on  board  a  Skylark  rocket  flown  from 
Woomera,  Australia  (31°S,  137°E,  magnetic  dip  63°)  at  1351  hours  LMT  un 
September  20,  1962.  Hall  and  Fooks  have  taken  into  account  Molmud's  [1959] 
correction  for  the  energy  dependence  of  the  collision  frequency.  At  88  km  the 
collision  frequency  was  calculated  from  the  ratio  of  the  amplitudes  of  the  up-going 
and  reflected  parts  of  the  wave  transmitted  from  a  ground  based  transmitter.  A 
curve  was  then  drawn  through  this  point  at  88  km  and  made  proportional  to  air 
pressure  taken  from  CIRA  1961 .  Hall's  and  Fook's  values  of  Hu  must  be 
multiplied  by  2/3  to  give  l/m  i.e.  ,  the  collision  frequency  for  mono-energetic 
electrons  of  energy  KT.  The  values  of  V m  so  determined  are  given  in  Fig.  9. 

Other  pertinent  data  not  shown  here  are  those  of  Saiah  and  Bowhill  [’966], 
Kane  [1962],  and  Landmark  [1965]. 

We  choose  the  data  of  Hall  and  Fooks  [1965]  for  our  model  ionosphere.  These 
values  are  given  in  Table  5. 

4.2  Electron  Cross  Sections  for  Momentum  Transfer  with  Neutral  Gases 

Inaccurate  results  yielded  by  theoretical  studies  of  low  energy  Mcctron- 
neutral  collision  cross  sections  have  necessitated  reliance  upon  experimental 
measurements  of  these  quantifies.  Banks  [I96bh]  has  made  a  review  of  such 
existing  data  and  has  chosen  approximate  expressions  for  these  cross  sections  for 
the  energy  range  appropriate  to  ionospheric  studies.  Hi  .  i  »sults  are  reproduced 
in  Fig.  10.  These  results  were  obtained  by  usmq  an  averaging  technique 
described  by  Banks  [1966b] . 

The  largest  cross  section  is  associated  with  H  and  that  of  N?  is  the  next 
largest  for  temperatures  greater  than  about  500°K.  The  cross  sections  for  He 
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and  O  are  essentially  independent  of  temperature  for  the  temperature  range 
of  interest.  Banks  quotes  the  uncertainty  in  Sen  (the  velocity  averaged 
electron-neutral  cross  section)  for  N2  and  O2  as  being  less  than  20  par  cent. 
Estimates  for  the  uncertainty  in  the  cross  sections  for  H  and  O  are  given  by 
Barks  as  +25  per  cent  and  +30  per  cent,  respectively.  The  cross  sections 
for  0^,  and  He  were  based  on  experimental  results  whereas  the  values 
for  O  and  H  were  derived  from  theoretical  calculations  of  scattering  phase 
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5 .  Concluding  Remarks 

No  attempt  has  been  made  in  this  report  to  suggest  ionospheric  production 
mechanisms  nor  to  deal  with  irregularities  in  ionospheric  structure.,  such  as 
sporadic  E.  We  have  also  ignored  the  problem  of  seasonal  variations.  Certainly 
there  are  seasonal  variations  in  electron  and  ion  densities  and  in  the  electron 
and  ion  temperatures.  There  may  also  be  some  seasonal  variations  in  the 
electron  collision  frequency,  although  this  effect  has  not  yet  been  resolved. 

Such  simplifying  assumptions  must  be  realized  when  using  the  data  contained 
herein  for  problems  having  a  seasonal  dependence. 

We  have  attempted  to  synthesize  existing  data  into  a  consistent  model. 

Where  the  experimental  data  are  lacking  we  have  relied  on  theoretical  predictions. 
Tor  the  electron  and  ion  density  profiles  and  the  electron  collision  frequencies, 
we  have  referred  primarily  to  the  data.  For  the  electron  and  ion  temperatures 
at  altitudes  above  1000  km  and  below  200  km  we  have  been  forced  to  emplc  . 
theoretical  values.  Also,  the  values  for  the  electron  cross  sections  for 
momentum  transfer  with  atomic  oxygen  and  atomic  hydrogen  were  based  on 
theoretical  results. 

Ion  collision  frequencies  and  cross  sections,  although  important,  have  not 
been  dealt  with  in  this  report.  These  parameters  are  considered  elsewhere 
[Shaeffer  and  Inoue  ,  1968b]. 
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On  the  following  pages  will  be  found  in  table  form  the  values  of  the  various 
model  ionosphere  parameters  obtained  from  the  adopted  curves  discussed  In  the 
previous  sections.  In  addition,  the  values  for  the  neutral  gas  densities  and 
neutral  gas  temperatures  from  CIRA  1965  Model  I  atmosphere  are  tabulated. 
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CHARGED  PARTICLE  DENSITIES  IN  THE  D  AND  E  REGIONS 
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TABLE  3. 

NEUTRAL  PARTICLE  DENSITIES 


ALT 

N(N’2) 

N(02) 

MO) 

N(Hc) 

MA) 

ALT 

N(N2) 

N(°2> 

N(O) 

N(HC) 

N(A) 

Knt 

cm -3 

cm"3 

cm-3 

cm“3 

cm”3 

Km 

cm" 

3 

cm"3 

cm"3 

cm"^ 

16 

3.60SE 

4 

1.606E 

14 
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2.950 

6.693 

3.210 
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51 

1.0(53 
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1.274 
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15 
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1 . 938 
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11 

1.757 

1.347 

2.934 
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62 

3.952 

1,060 

4.735 
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TABLE  4* 

ELECTRON,  ION,  AND  NEUTRAL  GAS  TEMPERATURES 


ALT 

Km 

To 

°k 

Ti 

°k 

Tn 

°k 
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Km 

Tc 

°fc 

Ti 

°k 

Tn 

°k 

52 
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265.5 
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54 
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TABLE  5 . 

ELECTRON  COLLISION  FREQUENCIES  AND  AVERAGE  CROSS  SECTIONS 

FOR  MOMENTUM  TRANSFER 
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94 
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96 
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3.25 
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98 
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U,  ABSTRACT 

-Full  vrave  solutions  have  been  found  for  plane  wave  and  localized  disturbances  in 
the  ULF  (micropulsation)  and  ELF  frequency  ranges.  Slow  and  fast  wave  disturbances 

incident  at  5000  km.  altitude  at  the  frequencies  100,  3.0,  1/3?  and  17300  ops.  were 
analyzed.  A  model  ionosphere  developed  from  existing  rocket  and  satellite  data  for 
venditions  of  midday  and  midlatitude  at  sunspot  minimum  was  used  in  t.he  calculations. 

The  -old  plasma  dispersion  relation  was  written  conveniently  in  terrs  of  the  localised 
difurlance  factor  and  the  ionospheric  conductivity  tensor  elements,  the  localised 
disturbances  were  found  to  produce  large  vertical  electric  fields  in  the  neutral 
atmosphere  for  the  incident  slow  waves  in  the  micropulsation  range  and  the  incident 
fast  wave  at  100  cps.  Localization  causes  strong  mode  coupling  in  the  lower  D  region 
a:;.i  the  coupling  coefficients  are  larger  than  the  quart ic  roots  at  1/3F0  cps.  Charge 
ox"».atige  of  C+  and  0  results  in  an  enhanced  Pedersen  conductivity  region  which  in  some 
ares  a  -i  s  as  a  cavity  wall  for  wave  energy  storage,  Charge  exchange  also  cause?  some 
i«,!  rorju  ion  at  Pci  frequencies  in  the  range  150-2'  0  ten.  The  conductivity  of  the  P.ar'h’s 
;• e  is  included  and  found  to  result  in  almost  '  otal  reflection  at  ground  level  for 
m ; '  frequencies  considered.  Incident  slow  mode  signals  at  3.0  cps.  undergo  polarize* ioi 
several  at  1000  km.  and  then  propagate  mostly  in  the  R  tvde  through  the  ionosphere. 

■  !»CC0  plasma  drift  velocities  between  0.1  and  ".0  tan/ser.  are  sufficient  for  the 
' idem  slow  waves  to  excite  typical  Pci  and  Pc5  ground  observed  signals. 
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